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Abstract 



Open Charm Production And Its Lepton Pair Signal 
In High Energy Nuclear Collisions 

Ziwei Lin 



Open charm production at the initial and the pre-equihbrium 
stages in nuclear collisions is studied. The pre-equilibrium con- 
tribution to open charm production in A + A collisions at ^/s — 
200AGeV is shown to be very sensitive to correlations between 
momenta and space-time coordinates of the minijets produced 
after the initial collision. A minimal correlation model based 
on the uncertainty principle is introduced, and leads only to a 
small pre-equilibrium charm yield as compared with initial yields. 
Opposite-sign lepton pairs ee,e/i and /i/i from open charm decay 
in p -|- ^4 reactions at ^/s = 200 AGeV arc proposed as a mea- 
sure of nuclear shadowing effects. Via an approximate scaling, 
the ratio of the spectra from p + A collisions to those from pp 
collisions is shown to reflect the shadowing function well. Finally 
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it is shown that the required measurements are feasible at RHIC. 
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Chapter 1 



Introduction 



The quark-gluon plasma phase of matter |jl|, |^ , a state consisting of weakly 
interacting deconfined quarks and gluons, was predicted to exist soon after 
the discovery ^ of the asymptotic freedom property of Quantum Chromody- 
namics (QCD). That novel phase of matter was predicted to occur perhaps 
deep in the cores of neutron stars and during the first few moments (mi- 
croseconds) of the Big Bang [|5[. It was also suggested that perhaps the 
quark-gluon plasma could be formed in the laboratory in very high energy 
collisions of hadrons or nuclei 0. Since the mid 1970's the experimental 
search for new phases of nuclear matter |^ using heavy ion reactions has 
been underway. Until 1984 only light ion reactions were possible to study 
at energies below 2AGeV (GeV per incident baryon) p. In the past ten 
years, experiments with heavy ion with energies up to 200 AGeV have been 
investigated [1^, |Tl], [1^. In 1999 a new era of experiments with heavy ion 



beams in the collider mode at the Relativistic Heavy Ion Collider (RHIC) at 
the Brookhaven National Laboratory (BNL) with center of mass energies up 
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to lOOAGeV will begin. 

In this chapter, we will review key elements of QCD thermodynamics and 
the dynamics of heavy ion reactions in order to put the thesis work described 
in the later chapters into context. The work on open charm that forms this 
thesis is related to the search for the QGP state of matter at RHIC energies. 
It deals with a probe of the earliest moments of heavy ion reactions and is 
expected to be one of the key diagnostic tools of that plasma. In Chapter [TT^ 
we will review, from both the perturbative QCD and lattice QCD calcula- 
tions, the basic idea about the phase transition from hadronic matter to the 
quark-gluon plasma. In Chapter we will introduce nuclear collisions as 
a tool of studying the quark-gluon plasma, and estimate thermalization and 
chemical equilibration processes in heavy ion collisions. Then in Chapter |1.3| 
we will review several of the important proposed observables including Heavy 
Quarkonia, dileptons, direct photons, disoriented chiral condensates (DCC), 
transverse flow, and Hanbury-Brown Twiss (HBT). Finally in Chapter |1.4| 
we will give the introduction on open charm signal, which is the main subject 
of this thesis. 

1.1 Quark- Gluon Plasma 

One of the most remarkable properties of the Standard Model of strong 
interactions, QCD, is that the fundamental color degrees of freedom carried 
by quarks and gluons can not be directly observed. That property, called 
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confinement, restricts physical observables to composite color neutral objects 
called hadrons and nuclei. It is thought that this property is caused by 
the very complex, non-perturbative nature of the physical vacuum. That 
vacuum structure is also thought to account for the observed breaking of 
the dynamical chiral symmetry of QCD The naive perturbative QCD 



vacuum is not the ground state, but has an energy density B ^ 200GeV/fm^ 
above the physical vacuum that is filled with gluon and quark condensates. 
The rich hadron spectroscopy and complex strong interaction phenomenology 
is caused by that non-trivial vacuum structure. As T.D. Lee has repeatedly 
emphasized over the past 20 years |T^, understanding that structure is one 
of the fundamental problems of modern physics. 

One of the few untested and key predictions of the Standard Model is 
that the structure of the physical vacuum could drastically change at high 
temperature and baryon density. Due to asymptotic freedom, as discussed 
below, the condensates could melt at sufficiently high temperature, and thus 
hadronic or nuclear matter is expected to transform into a weakly interacting 
gas of partons. 

1.1.1 Asymptotic Freedom 

From perturbative QCD it is known that the coupling constant depends on 
the momentum scale Q [Q. For SU{N) gauge theory with Uf flavours of 
fermions, this is expressed as the renormalization group equation: 

-^ = boal + hal + 0{at) , (1.1) 
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where t = ln((5^/A^), and 



^0 - ^-^f^-^) , (1.2) 



bi 



3AN^ f im 



+ — "/ 



1.3) 



For example, for QCD SU (3) theory the lowest order expression of the cou- 
pling constant is 

where Q is the momentum transfer and A ~ 200MeV. 

Therefore QCD has asymptotic freedom in that for < 6 the running 
coupling constant approaches zero at the limit oi Q —>■ oo. This raises the 
possibility that we will have a free gas of quark and gluon plasma when the 
nuclear matter is at high temperature or at high density 

1.1.2 Deconfinement and Chiral Restoration 

From the asymptotic freedom of QCD one expects to have an ideal gas of 
quarks and gluons at high temperature. However at low temperature quarks 
and gluons are subject to color confinement, and chiral symmetry is spon- 
taneously broken. Therefore at some high temperature we expect to have 
deconfinement and chiral restoration. 

The QCD Lagrangian can be written as 

£ = + J2i:,{D + m,) , (1.5) 
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where following quantities expressed in terms of the gauge field are: 

= d^Al-d^A^^ + gr'^A^Al , (1.6) 

D = i,(td, + ^gA';,X'^^ , (1.7) 

A, ^ IaI^X" , (1.8) 

where A° are the generators of the SU (N) group in the fundamental repre- 
sentation. The A*^ satisfy the following relations with /"'"^ being the antisym- 
metric SU{N) structure constants: 



2ir''X' , (1.9) 



Tr (A'^A'') = 25"^ . (1.10) 

The above Lagrangian in eq.( p..5|) has exact symmetry only when the 
quark masses are infinity or zero. In the former case we have the pure gauge 
theory with Z{N) symmetry, and in the latter case we have chiral symmetry. 

When the quarks are infinitely heavy, we can study the deconfinement 
phase transition from the breaking of the global Z{N) symmetry in the pure 
gauge theory. The order parameter is found to be the Polyakov loop operator 

L{x) =TrVexp(^i dtAo{x,t)^ , (1.11) 

which takes the trace of the path integral of an exponential function. 

The thermal expectation value of the Polyakov loop operator is related 
to the free energy of inserting a static quark in the vacuum: 

<L(f)>=exp(-^^) , (1.12) 
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where Fq stands for the free energy of the vacuum. 

In the confining phase the field can not be screened in the pure gauge 
theory, and the string potential grows linearly with distance, it therefore 
takes infinite free energy to insert a quark and < L{x) >= 0. However in 
the deconfining phase the field can be Debye screened and this results in 
< L{x) >7^ 0. Thus < L{x) > is the order parameter for the deconfinement 
phase transition and provides information on the breaking of Z{N) symmetry 
in pure gauge SU{N) theory. 

It was pointed out by Svetitsky and Yaffe [0 that the four dimensional 



SU{N) pure gauge theory at finite temperature is in the same universality 
class as the three dimensional Z{N) spin theory. Numerical studies on the 
effective models give strong evidence that SU{3) pure gauge theory goes 
through first order deconfinement phase transition [ll6| . 

The situation with dynamical quarks is quite different. The Z{N) global 
symmetry is explicitly broken by the finite quark mass. In general it is ex- 
pected that with dynamical quarks the confinement phase transition weakens, 
since the field can be screened by qq pair productions and the strings can be 



broken even in the confining phase [17 



On the other hand, when the quarks are massless, the QCD Lagrangian 
is invariant under the following transformations on the fermion fields ipf. 

^^ ^ e'"i^, , (1.13) 
V'. ^ e^"^^^, , (1.14) 



R 



;i.i5) 



where 75 is the Dirac matrix, and 



U^,U^ E SU{n 



;i.i6) 



1 + 75 
2 

1-75 



;i.i7) 
;i.i8) 
;i.i9) 



Hence the Lagrangian has the following chiral symmetry 



Uv{l) X Ua{1) X SUR{nf) x SUL{nf) . 



;i.2o) 



However, Ua{^) symmetry is broken by quantum effects, and the correspond- 
ing axial current [T^ 

4 = E^W5^. (1.21) 



is not conserved. Instead we have the Adler-Bell-Jackiw anomaly []T9 



as 



n 



f9' 



As a result the chiral symmetry is EC 



1.22) 



ZA{nf) X SURinj) x SUL{nf) 



(1.23) 



Pisarski and Wilczek |?T] studied the following Lagrangian with the above 
chiral symmetry 

Cch = -Tr (a^^^a''^) - -m|Tr ($^<l>) + c (L>et$ + L>et<l>^ 



Tr f^^^ 



2 VT^ 



(1.24) 



where the fields are 



^ij (x< 7Pi{l+-f5)^j > , (1.25) 

and where gi, g2 are related to /52 in the renormalization group. 

The order of the chiral phase transition from the above general Lagrangian 
is found to depend strongly on n/. The transition is first order for n/ > 3. 
For Uf = 1 the transition does not even exist. For the interesting Uf = 2 
case, the transition depends sensitively on the parameter c, the magnitude 
of the mass term. The phase transition is first order if c = 0, weakens when 
c becomes bigger, and even disappears if c is too large. The exact nature 
of the phase transition has also been studied with mean field analysis and 
lattice Monte Carlo calculations. 



1.1.3 QGP Perturbative Equation of State 

For QCD with massless quarks, the thermodynamic potential was calculated 
as a function of temperature 



-1 



InZ 



ao + 02^^ + a^g^ + (04/ In^^ + 04) g'^ 
+ (as; In/ + 05)/] + 0(/) , 



;i.26) 
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where Z is the partition function. The first four coefficients in the above 



equation are derived in terms of number of colors and flavours p2 
ao 



TT^ / 7 
— \N'^-\\ -Nn 
45 V 4 



02 



as 



144 



l(iV=-l) (iV + 



4 



97r V 



127r 



3/2 



(1.27) 
(1.28) 

(1.29) 
(1.30) 



327r^ V y V 3 6 
Then other thermodynamic quantities such as energy density, pressure 
and entropy density are expressed in terms of the thermodynamic potential 



P 



Tdv In Z 



T T dT ' 



(1.31) 
(1.32) 
(1.33) 



We expect there to be a free gas of quarks and gluons at the weak coupling 
limit (7 0, so the leading term oq in eq. ( |1.26| ) just gives the Stefan- 
Boltzmann quantities for massless bosons and fermions: 



TT 



P 



=90 



2{N^ -l) + -N ■2-2nf 
8 



3aoT^ = 3p , 



1.34) 



1.35) 



4anT3 



;i.36) 
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When one considers finite baryon density, one introduces the quark chem- 
ical potential /z, and then the equation of state (EOS) is given by [^^ 



P = 

^ 90 



2(iV2 - 1) + Ia^ . 2 ■ 2nj 
8 
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e = 3p , (1.38) 
n = ^(/i^ + vr^TV) , (1.39) 

s - '-^^ , (1.40) 

where n is quark number density. We mostly concentrate our discussions 
on the physics for zero baryon density case, since small baryon density is 
expected for very high energy heavy ion collisions. 

Perturbative methods provide us with a powerful tool to study QCD |^ . 



However, the result is an expansion in the coupling which is of the order 
1 at any practical energies despite the fact the it is a decreasing function 
of energy. Therefore the higher order terms in eq. (|1 . 26|) are comparable 
with the lower terms. Furthermore the expansion often encounters technical 
difficulties such as the summation of infinite terms with the same order. 

We take eq.( |1.26|) as an example. In order to eliminate the infra-red 
divergence at 0{g^), one has to introduce a magnetic mass in the gluon 
propagator p5| 



mM oc g'T . (1.41) 
Then there are infinite number of diagrams contributing at 0{g^), since a dia- 
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gram with n number of four gluon vertices gives a term which is proportional 

to 



^ • (1.42) 

mM J 



Ahhough eq.( |1.26|) is a perturbative resuh, it includes infinite resumma- 
tion of ring diagrams that take into account the polarization and screening 
effects of the plasma. Similar problems occur at all higher orders, and it 
seems impossible to solve the problem, which is called "Linde's disease" , an- 
alytically via perturbative expansion. Thus even when g is small, it seems 
that we still can not solve the problem analytically. Non-perturbative calcu- 
lations starting from first principle are therefore urgently needed. 

1.1.4 QGP Non-perturbative Equation of State 

In order to overcome the difficulties encountered in perturbative QCD cal- 
culations, lattice QCD was proposed by Wilson in 1974 pB|. Lattice QCD 



regularizes the QCD Lagrangian by fixing the lattice spacing. Then thermo- 
dynamic properties of the QCD can be studied by Monte Carlo simulation 
|27| . First one needs to discretize the QCD action 



fl/T r 

S= dxo C , (1.43) 
JO Jv 

where the QCD Lagrangian C is given by eq.( |1.5| ). 

Wilson proposed a way to map the gluon part of the action to lattice |2B 
First a link variable is defined to be 



Unn = exp 



;i.44) 



12 

where a is the lattice spacing, is a point in the four-dimensional Euclidean 
lattice, and /x labels the direction of the link. Then the gluon part of the 
action is 



where Re takes the real part of a function, and 



9^ 



1.45) 



(1.46) 



The above gluon action is equivalent to that in the QCD Lagrangian 
eq.( |1.5|) at the continuum limit a — 0. We expand the trace term in 
eq.( [OD to 0{a^): 



-iga{A^+A^)-iga'^d^A^+{igaf[A^,A^]/2 



;i.47) 



^ n+u,fi^ n,u — ^ 



igaiAfj, +A^)+igo? d„A^ + (iga)^ [A^,A^]/2 



'1.481 



TI TJ Tj'^ TI^ 



-iga2(0^A^-9^Ap-i9[A^,A^])+O(a3) 



(1.49) 



+ O(a^) 



G 



^JLy P 



1.50) 
1.51) 
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where 

V^ = {Nta){Nsaf = ^1"" (1.52) 

is the volume of the four dimensional lattice. 

The lattice regularization of fermion action turns out to have fundamental 
difficulties. The naive regularization produces 16 poles in the fermion prop- 
agator, thus results in 16 fermions from the every one original fermion. This 



fermion doubling problem can be solved in at least two ways. Wilson p6 
proposed to add a large mass term to eliminate all the other 15 unwanted 
fermions. However, the mass term explicitly breaks the chiral symmetry, 
which is very important for the study of QCD phase transitions. Kogut 
and Susskind proposed a way to reduce the 16 fermions to 4 (staggered 
fermion formulation), and this method has the advantage that part of the chi- 
ral symmetry still remains. Wilson's fermion method and staggered fermion 
method will agree and both produce the continuum fermionic Lagrangian in 
the limit a — 0. 

In the staggered fermion formulation, the fermion part of the action is 
given in terms of the fermion fields x{x)i which are anticommuting Grass- 
mann variables 

Sf = ^Xn,iQn,mXm,i ■ (1.53) 
i 

Here the fermion matrix Q includes the usual derivative term and the mass 
term. 
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Once one fixes the formulation of the fermion action Sp, the partition 
function on lattice becomes 

Z= fl[dUn,>.Y[dXn.,^dXm,^e-^''-''^ , (1-54) 

n,M m,i 

where Xm,i are the fermion fields on lattice. Then the expectation value of 
an operator O is 

< O >= ^l[dUn,^l[dXm,^dXm,^^e-^^-''^ . (1.55) 

Since the fermion action Sp is bilinear in the fermion fields, it can be inte- 
grated out and only the gluonic degrees of freedom remain: 

<d>=Z-^ III dUn,^d n det Q'e-^^ . (1.56) 



The lattice QCD results on phase transitions are summarized in Figure [O 
from ref.|^U[. The order of the phase transition is plotted as a function of 
the light quark mass m„ and the strange quark mass nig measured in units 
of a^^ . We proceed anti-clockwise from the upper- right corner. In pure 
gauge theory where all quarks are infinitely heavy, there is a first order phase 



deconfinement transition ||30[. With two flavours of quarks with small but 



finite masses, one has the two solid boxes on the upper border of Figure 1.1 



and there is no phase transition |2^, but there is a rapid crossover [^. For 
two massless quarks there is a second order transition. When we have three 
(or more) massless quarks, the transition is known to be first order; and 
it remains first order even for the case of three light quarks [^, which is 
represented by the lower solid circle. For QCD with one flavour of massless 
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quark, there is no phase transition. In lattice QCD physical quantities are 
expressed in terms of the lattice spacing a. To have a feeling of the magnitude 



of a, m^a ~ 0.4 in calculations of ref. I^Uj. Therefore the real world where 
one has two nearly massless quarks and one heavier strange quark lies in 
the area represented by the virtual circle in Figure |LT| . To probe near the 
real QCD, the lattice calculation with rriu^ci = 0.025, m^a = 0.1 (solid box 
in the middle of the figure) was performed but showed no phase transition. 
However, the current lattice calculations have uncertainties due to finite size 



effects. For example, in ref. |29| mj^/mp = 0.46 instead of the physical value 
of 0.64. Much more work needs to be done to determine the exact QCD 
phase diagram. 

The lattice QCD results on the equation of state of SU{3) pure gauge 



theory is shown in Figure |1.2| p2| . The upper dashed straight line is the 



Stefan-Boltzmann value for the energy density according to eq. (|1.35| ), and 



the other plotted physical quantities are the values extrapolated to the con- 
tinuum limit. One notices that the energy density quickly approaches the 
ideal gas value after T^, however the pressure approaches slowly and thus 
gives the biggest deviation from the ideal gas equation of state, e = 3p, at 
T ~ (1 — 2)Tc. From the perturbative QCD result given by eq. ( |1.26D , the 
deviation from the ideal gas equation of state is 

e-3p dn dg ^ 
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Figure 1.1: Phase transition from lattice QCD at finite temperature as a 
function of mu,da and m^a from ref. p9| . 
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where 

N(N^ - 1) 

Cl = ^-J^ C2 > . 

The deviation can also be connected to the gluon condensate 
e - 3p =< >o - < >T , 

where 

< >o~ 2 GeV/fm^ . (1.60) 

The zero temperature gluon condensate is important for the deviation of the 
quantity (e — 3p)/T^ when close to Tc, and its effect decreases rapidly at 
high temperature. The deviation at high temperature T/Tc > 2 could be 
attributed to the leading order perturbative result in eq. ( p,.57| ) assuming 
a large running coupling constant g'^{T) ~ 2. However the perturbative 
expansion can not be expected to provide us a quantitative result in this 
case. 

As seen above, the lattice calculation gives us valuable insight into the 
non-perturbative aspects of QCD. The critical temperature for the QCD 
phase transition appears to be about ISOMeV ||31|, which means a critical 
energy density of about IGeV/fm^ according to eq. ( |1.35| ) for three flavour 
QGP. However, one needs higher temperature (and thus a much higher energy 
density) for the onset of the perturbative behaviour, as can be seen from 
Figure |]2[ 



(1.58) 



(1.59) 
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Nevertheless, one has to be careful when drawing conclusions from lattice 
QCD results before one understands uncertainties such as finite size effects. 
To give an example, we show the results of the equation of state for two 
flavour QCD. In Figure |L3| [0^ energy density and pressure are plotted for 



two flavour QCD. The diamonds are the results on A^^t = 4 lattices, while the 
circles and squares are those on A^^^ = 6 lattices. The bottom straight line 
indicates the Stefan-Boltzmann value for energy density, and the other two 
straight lines indicate the Stefan-Boltzmann values corrected for flnite size 
effects. 



one 



By comparing A^^ = 4 results with Nt = 6 results in Figure |T]3 
flnds large flnite size effects. This is not unexpected since one has signiflcant 
flnite size effect in the free theory on lattices already. For example, the 
flnite temporal size effect for the energy density of an ideal gluon gas can be 
expressed approximately by [^| 



The flnite temporal size effect can therefore be as large as 50% on A^j = 4 



lattices for pure gauge theory. Improved actions on both gluon |^ and 
fermion part could help to limit the sensitivity to the flnite size. 

In order to improve the results from lattice QCD, larger machines are 
deflnitely needed. At present a 0.4 Teraflops supercomputer is under con- 
struction by N. Christ et al. at Columbia University |3^. That machine will 
be optimized for full QCD calculations, and will be one of the most promising 
lattice QCD projects to approach the continuum limit. 
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Figure 1.3: Equation of state for two flavour QCD on A^^^ = 4 and Nt = 6 
lattices as a function of T from ref . . 
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1.2 Why Do We Need Heavy Ion Collisions? 

In order to experimentally study the phase transition of the quark-gluon 
plasma, we need to achieve very high energy density in the lab. It is also 
important to have a large interaction volume in order to approach the ther- 
modynamic limit of the new phase. It is hoped that high energy heavy ion 
collisions can provide us with these necessary conditions. Since the initial 
interaction volume and parton multiplicity turn out to be large, the dense 
parton system may go through both thermalization and chemical equilibra- 
tion processes after its formation. At extremely high energies one may study 
an almost fully-equilibrated plasma of quarks and gluons with zero baryon 
density, and at currently accessible energies matter with high baryon density 
may be studies as well. The QCD phase transition is also expected to occur 
at high baryon density pf. 



1.2.1 Glauber Geometry 

In nuclear collisions many pp collisions occur almost simultaneously in a 
small nuclear volume. For hard perturbative QCD processes with small cross 
sections, the nuclear geometry and the QCD dynamics factorize to lowest 
order. Then it is assumed that they are independent two body reactions. 
For A + B collisions at a fixed impact parameter b, the number of binary 
collisions is given by ||37[ 



N^^{b) = J dpldyidy2d%d%S^'^\b - h - h 
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X Xl^„/A(xl,g^6^)x2^6/B(x2,Q^62)-^ . (1.62) 

For simplicity we do not consider nuclear shadowing effects in this section. 
Therefore the nuclear parton density in this case is simply 

^a/A{x, Q\ ^ = TA{^JalN{x. Q') , (1-63) 

where fa/N{x,Q^) is the structure function of parton a in the nucleon, and 
Ta{s) is the nuclear thickness function 

- /" dznA{y/z^ + \sf) , [ d^n^Ai^ = A . (1.64) 

J —oo J 

To show A scaling of hard processes, we take the simplest hard sphere 
model with radius Ra for the nuclear density function 

Ra = roA^/^ ro ~ 1.2 fm . (1.65) 

Therefore 



\s 



|2 



1-^ , (1.66) 



2nR\\ R\ ' 
^Ar^S(6) = j d%d%5^''\h-h-h2)TA{W)TB{h2) 

X j dp]_dyidy2Xifa/N{xi,Q^)x2fb/N{x2,Q^)-^{l.67) 

= TAB{b)aPP , (1.68) 

where 

TAB{b) = j d%d%5^''\h-h^-h2)TA{h)TB{h2) (1.69) 

= j d^^A{s)TB{b-^ . (1.70) 
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The scaling with atomic number A in central pA and AA collisions is, for 
example, given by 

r 3/4^/3 
Tp^(O) = J d'^,is}TA{s) ^ ^ , (1.71) 

^NP^/NPP oc A^/^ . (1.72) 

TMO) = / ci^^TK.-) = — , (1.73) 

^AT^^/ATPP oc . (1.74) 

An estimate shows that minijet gluons may be produced with a cross 
section ~ lOmb in pp collisions at ^/s = 200AGeV This estimate was 
obtained with a conservative transverse momentum cutoff Pq = 2GeV for the 
gluon scattering process. Thus, according to the above scaling with atomic 
number A, the number of minijet gluons may be on the order of thousands 
at RHIC. 

1.2.2 Global Features of Heavy Ion Collisions 

In this section we review the standard calculations for certain global features 
of heavy ion collisions such as energy density, temperature and multiplicity. 

At present there are two major facilities for heavy ion experiments: the 
Alternating Gradient Synchrotron at Brookhaven (AGS), and the Super Pro- 
ton Synchrotron at CERN (SPS). They are both fixed target machines. The 
AGS accelerates nuclear beams from proton to gold up to the energy of 
29{Z/A)GeV, and the SPS accelerates beams from proton to lead up to the 
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energy of 400(Z/A)GeV. Two larger colliders are under construction: the 
Relativistic Heavy Ion Collider at Brookhaven (RHIC) |^9||, and the Large 
Hadron Collider at CERN (LHC) [|0]. RHIC is scheduled to be completed 
in 1999 with the center of mass energy of 200AGeV for Au + Au collisions. 
Heavy ion physics is one important part for the LHC program, and the center 
of mass energy for Pb + Pb collisions at LHC will be 5.4ATeV. 

The total multiplicity rapidity density for pp collisions from SPS to Teva- 
tron energies has been fit with 

For AA collisions, that density can be parameterized approximately as 



,dy J \dy J 



where a can be determined at least from lower energy pA data 

a ~ 1.1 . (1.77) 

In the standard scenario based on Bjorken's longitudinal boost-invariant 
boundary conditions [^], the initial volume is assumed to be 

Vo = vrr^A^/Vo, tq ~ Ifm , (1.78) 

and is furthermore assumed to grow linearly with proper time. Given the 
average transverse momentum of final particles, 

< A;^ > ~ 400MeV , (1.79) 

^<m±> ~ SOOMeV . (1.80) 
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Hence the initial energy density is 

^^^ < ,„.>(.W- ^^^^..,3,„(^) «eV/to3 . ,,31) 

From this equation, we reahze that the initial energy density grows as a power 
of A, the atomic number of the heavy ion beam, but only grow logarithmically 
with the center of mass energy. Therefore it is much more efficient to reach 
a higher initial energy density by using heavy ion collisions than increasing 
the beam energy alone. 

According to eq. (|1.35|) we get the initial temperature assuming that the 



initial state is an ideal gas of gluons and three flavor of massless quarks: 

n=(^y^'^el^'xl50MeV , (1.82) 



,47.57r2 

where eo is expressed in the unit of GeV/fm^. 



In Figure 1.4 we plot the initial energy density and temperature as a 



function of the center of mass energy. We use eq. (|1.81| ) and eq. ( |1.82| ) from 



the standard scenario, and assume AA collisions with A ~ 200. 

Note that the multiplicity density dN/dy is about 300 times the energy 
density (in the unit of GeV/fm'^) according to eq. (|1.81|) . In Figure |L4 



we see 



that almost all four experimental facilities can achieve an initial energy den- 
sity above the critical ec — IGeV/fm^ for the QCD phase transition. However 
at the currently available AGS and SPS energies, the baryon stopping power 
is sufficiently large that matter is formed with high baryon density. We recall 
from the discussion of QCD thermodynamics that much less is known about 
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Figure 1.4: Initial energy density and temperature as a function of center 
of mass energy of heavy ion collisions. 
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the high baryon density case than the zero baryon density case, and we ex- 
pect to gain valuable knowledge from these experiments. At future collider 
facilities, matter is expected to be formed at much lower baryon densities. 

We have to stress that the above estimates for the global quantities in 
this section are based on a standard scenario, which only provides an order of 
magnitude estimate. Much higher energy densities may be achieved by taking 
multiple minijet production into account, as we discuss later in Chapter 



371]. Factors such as a bigger < k± > ||4^, a larger a from more secondary 



scatterings |44] as found in cascade models, and a smaller tq as in the hot 



glue scenario |^5[ can all result in higher temperatures for the initial QGP 
system. However, we note that even the conservative estimate P2[ stated 
above gives a high enough initial energy density to motivate the search for 
quark-gluon plasma in heavy ion experiments. 



1.2.3 Thermalization 

The quark-gluon plasma is defined to be an asymptotically free gas of quarks 
and gluons with both thermal and chemical equilibration. After the dense 
system of quarks and gluons is formed in high energy heavy ion collisions, it 
will approach thermalization via rescattering, and approach chemical equili- 
bration via creation and annihilation processes. However, the parton system 
only has a finite lifetime. Thus it is important to know the time scales that 
the system requires to reach thermalization and equilibration, and whether 
the parton system can achieve both stages in its finite lifetime. 
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In this section we first follow Shuryak's estimate of the above elastic 



scattering cross sections , and remark on the limitations of that approach, 
and then estimate the thermalization time. Finally we introduce new results 
on parton energy loss which imply that faster thermalization may be 
possible via inelastic processes. 

The cross sections for the relevant elastic processes are 



da 

lit 



M. 



99^99 



M. 



qZ 



9 / ut us st\ 
2V ~ ^ ~ ¥ ~ 11?) 

9 + + 



9q^9i 



qq^qq 



4 US 
4 + 

9 



t2 



(1.83) 
(1.84) 

:i.85) 

(1.86) 



One may estimate the above cross sections by separating them into large- 
angle and small-angle scattering parts ^5[. Since most entropy in high energy 
heavy ion collisions is in gluons P3], let us consider gluons. 

For an ideal gas of gluons at a temperature T, 



J 7r2 Jo e'-l 



2(iV^-l)C(3)^3^,^^3 



C(3) 



1.202 



2{N^ - l)7r- 
30 



T^ = a2T^ , 



;i.87) 
;i.88) 
;i.89) 
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= ~ 3T . (1.90) 

Ug 

For an ideal gas of quarks, equations become 

n, = 6lT^ b, = ^^4^ , (1.91) 



= hT\ &2 = ^ . (1.92) 



For large angle scatterings, the matrix elements in eq. (|1.86|) are finite. There- 



fore one might estimate the corresponding relaxation time using the average 
values of the Mandelstam variables for an ideal gas of gluons: 

<s>/2 = -<t>=-<M>~ {3Tf , (1.93) 

^ < M^^gg >~ 30.4 , (1.94) 

1 ^ 

ng < a >~ Ug <M^><t >~ 5a^T.(1.95) 



^large— angle 9 ^ < S >^ 



For small angle scatterings, the matrix elements are divergent att,u^ 0. If 
one just use the finite Debye mass |4^ to regulate the divergence [^, then 



the contribution from small angle scattering is 

tmin = f^l = g'T^ , (1-96) 

1 97ra^ 

- ~ 2.2a,T . (1.98) 



small— angle j. 
Tg l^min 

Therefore the estimate of the rate for gluon elastic scatterings is 

1 



~ 5a^T + 2.2a,T . (1.99) 



^9 
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The above method |^5| of treating the small angle divergence is far from 
exact. Although Debye screening eliminates the divergence from longitudinal 
exchange of gluons, the transverse interactions, i.e. the color-magnetic inter- 
actions, can not be simply screened. This problem was discussed earlier by 
Danielewicz and Gyulassy |^. Later Baym et al. ||4^ found that the Lan- 



dau damping of transverse gluons provides an effective cutoff at a distance 
~ for the infrared divergence of transverse interactions. The result for 
the gluon relaxation time considering the viscosity of a quark-gluon plasma 
in the weak-coupling limit is given by 

— ~ 4alln{l/a)T . (LlOO) 

^9 



Comparing eq. ( |1.99|) to eq. (|1.100|) , we note that the more precise depen- 



dence of the relaxation time Tg on the coupling constant a in eq. (|L100|) is 



hinted by the combination of linear and quadratic a dependent terms in the 
crude estimate found in eq. (|1.99|) . 

Now we take eq. ( |1.100D and find that for a high temperature T ~ O.SGeV, 
the thermalization time for the gluons is large: 

a ~ 0.3 , (LlOl) 

~ Ifm . (L102) 

Simila ranalysis gives a quark relaxation time which is about twice larger 
than the gluon relaxation time. Thus, via elastic scatterings it is difficult to 
quickly approach thermalization. 
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However, thermalization may be accelerated by inelastic processes such 
as parton radiations from multiple scatterings in the dense medium and/or 
gluon multiplication processes {2 g ^ n g, n > 3) |5^. The radiative en- 
ergy loss due to multiple scatterings was first considered in QED by Landau, 
Pomeranchuk and Migdal (the LPM effect) ||51|. Similar effects were con- 
sidered in QCD by Gyulassy and Wang and recently new calculations 
from Baier, Dokshitser, Mueller, Peigne and Schiff (BDMPS) yield un- 
expected results. Instead of being finite, the energy loss per unit length in 
QCD for a fast parton grows with the energy of the parton as ^/E when the 
medium is infinite, or grows with the medium length L when the parton is 
extremely energetic. The BDMPS result for the spectrum of radiative gluon 
energy loss per unit length is given by p6| 



where k = Ag/i^/cj, w is the frequency of the radiation, /i is the Debye mass 
induced by the medium, A^, is the mean free path for the gluon, and Cr is 
the Casimir operator for the parton. 

In terms of radiation frequency, the LPM effect in QCD occurs for 



while the Bethe-Heitler regime lies below ujbh, and the factorization regime 
lies above the frequency upact and up to the energy of the parton, E. 

In the case of an infinite medium, the LPM effect occurs for frequencies 




(1.103) 



UJbH ~ A/i < UJ < LUpact 




(1.104) 



A 
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up to the order of E, and therefore 



dE f ( dl \ f /I 

— = jduj (^^^j oc J dujJ-\nuj oc y/E\nE . (1.105) 



dz 



While in the case of an infinite energy E, the LPM effect stops at ujpact, 
thus 



dE 

— oc LlnL . (1.106) 

dz 



This nonhnear behaviour may imply a huge radiative energy loss for fast 
partons in a dense system. The strong radiation may accelerate the process 
of thermalization, which is probably slow when one only includes the elastic 
scatterings according to the previous estimate. 

The previous estimate of thermalization time, e.g. eq. ( |1.99| ), assumed 



at the initial stage a chemically equilibrated gluon gas. With a high ini- 
tial temperature, the corresponding gluon density is expected to be large. 



eq. ( 11.881) . However, the initial gluon density produced via minijets and soft 
mechanisms may still not be high enough to populate the whole phase space. 
In order to study the chemical equilibration of the dense parton system, as 
will be discussed in the next section, one needs the initial condition including 
thermalization time scale, temperature and chemical compositions for quarks 
and gluons. Therefore one has to take the non-zero chemical potential into 
account. 
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1.2.4 Chemical Equilibration 

Chemical equilibration is approached via flavour creation and annihilation 
processes. The evolution of chemical composition can be described by a set 



of master equations |£3|. In order to solve the equations one usually assume 
that the partons have reached approximate thermal distribution. As we will 
see, the time scale for chemical equilibration is much longer than that for 
thermalization, and therefore the above assumption is supported. 

As discussed in the previous section, gluons and quarks in general are not 
in chemical equilibrium. With non-zero chemical potential, one can introduce 
fugacity parameters and Xg 



/-EAre—'^^rS^-J.^. (1.107) 

where i = g,q, (3 = 1/T, and u'^ is the velocity of the local co-moving frame. 

We study the following dominant processes for parton creation and anni- 
hilation 

99^ 999, 99 ■ (1.108) 

For a baryon symmetric system where = n^, the evolutions of the parton 
densities are 



1 n Ug 



dAn,un = 2^2n^l-^l , (1.110) 
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where 

c^s = < a{gg ^ ggg)v >, a2 =< a {gg qq) v > , (l.Ul) 



n, 



9 



Ug/Xg, Ug = Ug/Xg . (lAU) 



The temperature evolution is 

d^{eu'')+pd^u'' = . (1.113) 

During early times one can neglect transverse expansion, and assume the 
Bjorken longitudinal expansion, hence the above equation becomes 

^ + i±i^ = 0. (1.114) 

dr T 



For a thermally equilibrated quark and gluon gas, according to eqs. (|1.90|) 
and ( |1.92| ), the energy density becomes 

e = 3p = {Xga2 + Xgb2)T^ . (1.115) 

Hence the solution of eq. ( |1.113| ) is just 



= (A^as + Agfea) TV^/^ = constant . (1.116) 
With the assumption of pure longitudinal expansion, one has 

d,in,un = ^ + ^ . (1.117) 

OT T 

Therefore one can write eqs. ( |1.109| ) and ( |1.110D as 



2^ 



^ + 3^ + ^ = (1 - A,) - 2i?2 (l - ^ j , (1.118) 
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where the overdot represents the derivative with respect to r, ai and bi are 
from eqs. ([TSSD and ( [Ol , and 



1 1 



1.120) 



The above rates -R3 and R2 contain infrared divergences. One can estimate 
the rates using the Debye screening mass for gluons and the thermal quark 
mass to regulate the divergences H3, ISH]. Numerical results for the esti- 



mate are shown in Figure. 1.5 |55 



The initial conditions for above differential equations as in the case of 



Au + Au collisions at RHIC energy ^/s = 200AGeV is 



Tiso = 0.7fm, no = 3.2fm-^ < >= 1.17GeV . 



Using eq. ( pT^ ), eo = 4no <k^> /n, and A° = 0.16A° |55 



;i.l21) 



n = {Xgtti + Xqbi) . 



One can obtain initial values of fugacities and temperature 



= 0.05, A^ = 0.008, To = 0.57GeV 



1.122) 



1.123) 



Given the above initial conditions, the chemical equilibration eqs. ( |1.116|) , 
( |1.118|) and (|1.119|) can be solved numerically, and results of which are shown 
in Figure |T]^ |^ . 



In Figure jO] the temperature in this calculation drops faster than in 
Bjorken's free streaming case, since some energy is used to produce more 
partons in order to approach chemical equilibration. Gluons are much more 




Figure 1.5: The gluon and quark production rates from ref. 
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Figure 1.6: Chemical equilibration in Au + Am collisions at RHIC from 



ref . \EE . 
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likely to approach the chemical equilibrium than quarks. However, both of 
them can hardly reach chemical equilibrium within the lifetime of the plasma 
state. 

One can see that the result in this study depends strongly on initial 
conditions. This particular set of initial conditions may be conservative. As 
discussed in Chapter 1.2.3| , if fast thermalization could be achieved via inelas- 



tic processes, then one would begin with high initial parton fugacities, and 
the lifetime of the plasma would become longer. All of these contributions 
would accelerate the chemical equilibration. However, initial parameters are, 
unfortunately, quite uncertain since we have insufficient understanding of 
the soft physics. There are uncertainties in cutoff scales, rescattering cross 



sections, 2g ng process [Q, and small-x gluon distributions. Therefore 
careful investigations are still needed to determine whether more equilibrated 
plasma can be formed at RHIC 

1.2.5 Phase Diagram 

From lattice QCD numerical simulations, the QCD phase transition happens 



at Tc = 150 ± lOMeV [^] at zero net baryon density. The QCD phase 
transition is also expected to occur at some high baryon density nf at zero 
temperature, when the hadrons are so dense that the partons from different 
nucleons start to overlap. At present, there is no lattice QCD calculation on 
phase transitions with finite baryon density, and therefore there is uncertainty 
with the critical baryon density . In general one expects that the QCD 
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phase transition occurs on a smooth curve on the plane of temperature T and 
the baryon chemical potential /i^- In terms of baryon number densities, there 
is a jump from the hadronic phase to the QGP phase, just like the jump in the 
energy densities at zero baryon density. The phase diagram for the transition 
from the hadronic phase to the QGP phase is roughly plotted in Figure [l.T] 
as a function of the temperature T and the baryon chemical potential /i^. 
Figure p..8| shows the phase diagram as a function of the temperature T and 
the scaled baryon number density ns/n^^ where ~ 0.14/fm^ is the normal 
density in nuclear matter. The shaded area in Figure |1.7| represents the error 
from our current knowledge about the phase transition, and the shaded area 
in Figure |1.8| mostly represents the mixed phase. 

Let us estimate the critical chemical potential at T = 0. For high density 

quark matter consisting of m, d quarks at zero temperature: 

12 /-M, „ 2 o , , 

(fp=-nl , (1.124) 



(27r)3 7o 7r2^'^ ' 

^ ■ (1-125) 



27r2^'? 

If we assume that the QCD phase transition at zero temperature needs ap- 
proximately the same critical energy density as in the zero baryon potential 
case, then 

ec ~ IGeV/fm^ , (1.126) 

/o 2 \ 1/4 

^11^= f^^J ~ 480MeV , (1.127) 
^ ~ 3/i' = 1.4GeV , (1.128) 
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Figure 1.7: Schematic QCD phase transition diagram as a function of T 
and /iB- 
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rf 



= 0.93/fm^ ~ 7n^ in quark matter phase. (1.129) 



These estimates are consistent with the zero temperature values in Fig- 
ure [TT^ and Figure |0[ In order to obtain a general estimate on the phase 
transition boundary, expressions for the pressures in the two phases need to 
be known. The pressure Pp{T,fi) in the plasma phase is derived perturba- 
tively and expressed in eq. ( |1.37|) . The hadronic equation of state can not 
be obtained from perturbative calculations, but can be parameterized for 
different hadron species. In this simple case one gets the parameterization of 
the EOS in the hadronic phase by adding the parameterizations for nucleons, 
pions and other resonances together. The energy density for nucleons, for 
example, can be parameterized as 

K{n/no - If 3T" 



n 



m -Wq + 



;i.i3o) 



18 2 

where n is the nucleon density, m is the nucleon rest mass. Wo is the binding 
energy, K represents the compressional energy density, and the last term is 



the thermal energy [^]. One then take Gibbs' criteria of phase equilibrium 
Ph{T,fi) = Pp{T,fi) to determine the phase boundary curve Tc{fi), which is 
schematically shown in Figure |1.7| |]57| . 
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1.3 Signals of the Quark-Gluon Plasma 



In this section we discuss several key observables of the quark-gluon plasma, 



as sketched in Figure [1.9| . Among the probes of the early formed quark- 
gluon plasma, electromagnetic probes such as direct photons ^ and 
dileptons |6^ are the most direct signals, since they do not experience the 



subsequent thermalization and final state interactions. Heavy quarkonia |61 



and heavy quark fi^ , |62| productions are also proposed to be the sensitive 
probes of the properties of the QGP. Recently new interesting experimental 



results such as low mass dilepton enhancement p3| , |64| and J/ip suppression 
m A + B collisions were presented. These results were heavily discussed 
at the Quark Matter'96 meeting ||T2|. 



1.3.1 Heavy Quarkonia 

The production and suppression of heavy quarkonia bound states, such as 
J/ip, was proposed by Matsui and Satz [^] to be an ideal test of the quark- 
gluon plasma production. The authors claimed that "there appears to be 
no mechanism for J / ip suppression in nuclear collisions except the formation 
of a deconfining plasma, and if such a plasma is produced, there seems to 
be no way to avoid J/ip suppression." The proposed suppression signature 
of QGP formation was first observed in 200 AGeV O + U reactions at the 



CERN/SPS in 1987 . However, far from settling the matter, the debate on 



the interpretation of those data and on the uniqueness of the QGP mechanism 




Figure 1.9: Several key observables of the quark-gluon plasma. 
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for suppression has only intensified since tlien. 

The nuclear suppression of charmonium was observed in p+A reaction and 
was found to depend linearly on A^^^. In addition, the suppression was also 
observed for ip'. Since the QGP production is noi expected in p+A collisions, 
alternate mechanisms for suppression of quarkonia production were proposed. 
One attempt to explain the p + A data was to ascribe the suppression to 
nuclear gluon shadowing . (We will discuss more fully the physics behind 
this effect later.) However, this explanation was later excluded because the 
lowest energy data at 150GeV did not fit into the universal shadowing curve. 
The most satisfactory explanation thus far |^8[| of J/ip suppressions in both 
p + A and light ion O + U and S + U reactions is that it is due to pre- 
formation inelastic final state interactions of the color-octet ccg component 
which is produced in the initial hard glue-glue fusion event, with nucleons 
in the target with a cross section ct^at ~ 6mb. In this nucleon absorption 
model, the J/ip survival probability in nucleus A can be estimated as 



Sa = exp 



^0 



exp{-a^NPoLA) , (1.131) 



where z is the longitudinal coordinate along the beam axis, and b is the 
impact parameter of the initial gg — > cc hard pQCD event in the target 
nucleus with density profile, pA, normalized to A. The survival probability 
is the probability of having no final state interactions in the remaining part 
of the nucleus. It depends exponentially on the effective path length. La, 
travelled by the cc pair in a nucleus of mean density po ~ 0.14/fm^. We 
note that typically the average rapidity of the produced cc pair is high in the 
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target frame and that is why the straight hne (Eikonal) approximation can 
be used. 

In A + B colhsions, the data on quarkonia production are usually sorted 
in terms of an external calorimetric transverse energy trigger, Et- The trans- 
verse energy is usually measured in a window with restricted rapidity and 
azimuthal angle coverages, and depends sensitively on the particular exper- 
imental device and configuration. Generally only a detailed Monte Carlo 
nuclear event generator such as HIJING |70| or VENUS [71] can provide a 



way to calculate the dependence of the observed transverse energy on the 
impact parameter of the collision. With such a generator, the conditional 
transverse energy probability density, P^^{Et', b), that an A + B reaction at 
impact parameter b will result in a measured Et, can be estimated. 

The survival probability of J/ip in an A + B reaction triggered with a 
given Et can then be calculated as 



a^yi^T) 



oc -3^ / d^P^^iEr, b) dzAd%A dzBd%B S\b-bA - bs) 
AB J 

X pA{ZA,bA)pB{ZB,bB)e ■''A ■'-B 

— ^-Ci>NPQLAB{ET) ^j^ 132) 



In the first line the Drell-Yan cross section [|72| is used to divide the quarko- 
nium cross section, since the absence of final state interactions in qq —>■ fi^fj^~ 
makes the atomic number scaling of that cross section to be particularly 
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simple. Integrated over transverse energy triggers and hence impact pa- 
rameter, — ABa^Y ■ Fo^' the most central collision (high Et) trigger 
N^Y ~ A'^/^N^Y up to nuclear shadowing corrections. We note that both 
the above quarkonia and Drell Yan cross sections are understood to refer to 
comparable rapidity and transverse momentum kinematic regions. 



In Figure |1.1CI| the ratio of J /if) to Drell- Yan cross sections is plotted 
against the effective nuclear absorption length L = {Lab{Et)) ■ Unfortu- 
nately, the L dependence on the experimental observable Et is rather model 



dependent . Nevertheless, the exponential fit shows that both p + A and 
S + U data can be well explained taking a^N ~ 6 mb. 

A caveat in the above nucleon absorption explanation of suppression in 
light ion reactions is that a larger suppression of ip' is observed m S + U colli- 
sions than for ip . This can only be understood in terms of an additional 
suppression mechanism that only applies to the larger sized ip' states in the 
co-moving high density hadronic matter formed along side the quarkonium 
state. 



The most provocative discrepancy seen in Fig. 1.10, however, is the failure 



of the nucleon absorption model to explain the latest data showing enhanced 



suppression of J j^' in NA50 Pb+Pb experiment ||6^. The observed enhanced 
extra suppression of J/ ip in Pb + Pb data has been interpreted as the 
first direct evidence for the formation of quark gluon plasma. 

The above claim has however been challenged by results of a recent study 
by Gavin and Vogt based on a more careful analysis of the co-mover sup- 
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Figure 1.10: J/tp to Drell Yan ratio in p + A, S + U and Pb + P6 collisions 
as a function of the average path length L from ref . |7^ . 



48 



pression model . In an earlier study of NA38 S + U data, it is argued 
||78| that nucleon absorption alone (even with a cross section of 7mb) can not 
explain the Et distribution of S + U data, while co-mover dissociation cross 
section of Smb together with nucleon absorption of 4.8m6 could. Using the 
same parameters and updating the kinematical variables appropriate to the 
NA50 experiment, this model predicted the extra suppression of in Pb + Pb 
in quantitative accordance with the recent data. 

The J/ip survival probability due to co- mover scatterings can be simply 
estimated according to 



Sco = exp(^-J dTn{T)acoVre}j , (1.133) 

where due to Bjorken longitudinal expansion, the local co-moving density 
decreases with proper time as 



n(r) = no— . (1.134) 

r 



Here, tq is the formation proper time of the quarkonium state in the frame 
where its longitudinal rapidity vanishes. The initial co-moving density is 
found to scale approximately linearly with transverse energy, tiq = uqEt / 



| 77| . The final proper time, tf, marking the time that the ip exits the in- 
teraction it propagates in the transverse direction is approximately 
given by 

TF^ — . (1.135) 

Vrel 
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The survival probability after the co-mover scattering with cross section (Jc 
introduces another exponentially decreasing factor with increasing Et'- 



Sco ^ exp 



-acoVreinoToln f-^) Et/Et 



exp{-pET) . (1.136) 



Therefore, together with the nucleon absorption factor the survival probabil- 
ity including co-mover dissociation processes is given schematically by 

Sab{Et) oc < SaSbSco >~ e-'^^^o^-^^^^^^e"^^^ . (1.137) 

Although the average nuclear absorption length Lab{Et) grows with Et, 
it has a geometrical upper limit ~ (Ra + Rb) in central collisions. Therefore, 
nucleon absorption mechanism saturates at high Ex- However in central 
collisions Et may fluctuate to higher and higher values as the local rapidity 
density of hadrons increases. Co-mover dissociation can therefore increase 
indefinitely with increasing E^- The result of the detailed analysis [ [T? 



IS 



shown in Figure |1.11| . The suppression of J/ip production as a function of 



the actual experimentally observed Et is plotted for NA38 S + U [|^, |74 



and NA50 Pb + Pb data [^. Furthermore, the co-mover model could also 
be adjusted to fit well the observed ip' data in p + A, S + U and Pb + Pb 
collisions. This agreement indicates no compelling need to invoke the QGP 
explanation at this time from purely the phenomenological point of view. 
Theoretically, the validity of the large co-moving dissociation cross section 
has been questioned [Q, but from the phenomenological side all the available 
data seem to be well accounted for by this conventional mechanism. 
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Figure 1.11: Co- mover J/ip suppression vs S + U and Pb + Pb data from 
ref. [Ifl. 
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Another ambiguity associated with the interpretation of quarkonium data 
in terms of QGP formation is due to uncertainties in the basic production 
mechanism of such bound heavy quark state even at the pp level. In the 
generalized version of the color evaporation model ||8l|, it is assumed that 



charmonium cross sections are proportional to the sub-threshold charm cross 
section. With a fit proportionality constant of 3% P2[, the model seems 



to give a good description (with a factor of 2) of the energy dependence of 
the total J/ip production with pion and proton beams. It is also consistent 
with ip' /{J/ip) ratio data. This model has only one parameter for any bound 
state. However, it remains to be justified that this factor is independent 
of processes. In particular, the factor for pion and proton beams differs 
somewhat and the color evaporation mechanism can easily depend on the 
environment and hence on A in nuclear reactions. 



In the color-singlet model p3|, on the other hand, it is assumed that 



charmonium cross sections are proportional to the cross section of forming 
a color-singlet cc state with the corresponding quantum numbers. Then the 
non-perturbative probability for the color-singlet state to convert to charmo- 
nium is related to the derivative of the wavefunction of the singlet state in the 
origin. This model was able to give absolutely normalized production cross 
section, and the ratio of different charmonium states. However, it completely 
fails to reproduce the high px J/ip data by orders of magnitude. Formally, 
its problem can be seen |S^,e.g. via the infra-red divergence of the P-wave 
charmonium cross sections. 
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Recently, it was proposed that in some region of the phase space color- 



octet components might be dominant Contrary to the color-singlet 



model, the short-distance part of the interaction could also involve a color- 
octet cc. Although the color-octet matrix elements are suppressed by powers 
of V, the relative velocity of the heavy quark in the quarkonia, color-octet 
mechanism might dominate in the processes where the leading term in v is 
suppressed by other small parameters such as strong coupling constant 
or mQ/p± at large transverse momentum. The conversion of the color-octet 
state to a color-singlet charmonium is described by non-perturbative matrix 
elements. With those color-octet matrix elements fixed by the data, the 
model gives good agreement with the shape of the p± distribution. Although 
the magnitudes of the unknown matrix elements are consistent with the non- 
relativistic QCD model (NRQCD) |S^, more theoretical work is needed. 

Given the above uncertainties in the production dynamics at the pp level 
as well as the ability for conventional nucleon absorption and co-mover dis- 
sociation mechanisms to explain the quarkonium data, it is important to 
consider a much wider class of observables in the search for QGP formation. 

1.3.2 Dileptons 

Dileptons are one of the direct electromagnetic probes. Thermal dileptons, 
for example, can tell us about the thermal profile of the dense plasma. In ||59 



it was proposed that in the hot-glue scenario, where the plasma is less chem- 
ically equilibrated but hotter than in the chemically equilibrated plasmas, 
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dileptons of a few GeV energy will be enhanced over the standard scenario 
by a factor of 2. 

However, in order to observe the interesting thermal signals, one has to 
deal with a large dilepton background. One such background is the dileptons 



decayed from heavy flavor meson pairs. In a study by Vogt et al. , it was 
suggested that dileptons from open charm decay would be about an order of 
magnitude above the thermal and Drell-Yan signals in the few GeV region, 
and therefore it would be very difficult to subtract this large background 
and observe the desired plasma signal. However, Shuryak ||8^ suggested that 
one important effect was missing in the arguments of that study [^], which 



was the energy loss of the charm quark in the dense medium formed after 
high energy A + A collisions. Once one assumes that the charm quark would 
lose a finite amount of energy along its way (typically 2GeV/fm), then the 
dilepton spectrum from open charm will be pushed to the low mass region, 
and dileptons with a few GeV energy would be suppressed by over an order 
of magnitude. More theoretical work on the energy loss in a finite nucleus 
needs to be done before we know the magnitude of the suppression on high 
mass dileptons in A + A collisions. 

Very interesting data on low mass dilepton enhancement was reported at 
the Quark Matter'95 meeting [p^l ODI. It was found that dileptons in the low 



mass region 0.2 — 1.5 GeV are enhanced by a factor of 5 in + Au collisions 
In Figure I.V2 the CERES S + Au data is presented [^|. Similar effects 



were also observed in S + W collisions [pT |. 
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Figure 1.12: Low mass dilepton enhancement from CERES experiment from 
ref. [Bl. 
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Dileptons from vector mesons are of interest |^2|, ^ because they can 
provide information on the possible variation of hadron masses due to chiral 
symmetry restoration. The low mass dilepton enhancement was also one of 



the focal points of the Quark Matter'96 meeting [12 



It was proposed [|^, |9^ that the CERES data |^ could be explained by a 



drop of p meson mass in nuclear medium. In Figure |1.13| the assumption |53 
of in-medium mass reductions of nucleons and p mesons is plotted for several 
baryon densities. The p mass was assumed to have dropped to 270MeV at 
the density 0.4/fm^ and temperature 165MeV. However, this assumption of 



p mass drop is too large, and is inconsistent with QCD sum rules [Q. 

In Figure |1.14| the CERES data on low mass dilepton invariant mass 



spectra is plotted vs calculations from ref . ||95[ assuming the in- medium meson 
masses assumed in Figure |1.13| . With the significant p mass reduction in 
the dense hadronic matter, the agreement with the CERES data is greatly 
improved in the invariant mass region from 2171^^ to nip. 



In calculations from another group ^ a smaller p mass dropping was 
assumed: 



m* ~m°(l-0.18pB/po) 



1.138) 



and these calculations also explained the CERES S + Au data, as well as the 
HELIOS S + data gT[. 

We see that the explanations of the observed low mass dilepton enhance- 
ment are not unique. Different scenarios of p meson mass shifts in medium 
were assumed. In addition, the conventional scenario of the broadening of 
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Figure 1.13: Assumption of in- medium nucleon and p meson masses as a 
function of temperature from ref. |^ . 
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Figure 1.14: CERES data on low mass dileptons vs the calculation from 
ref. [0 assuming in-medium meson masses. 
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the p meson spectral function in medium |^ could enhance the low mass 



dileptons. The present experimental error bars are still too large to pin down 
different models, and therefore more precise experimental data and further 
theoretical studies are needed to extract the physics from the CERES and 
HELIOS data. 

1.3.3 Direct Photons 

Direct photons are also one of the direct probes which could provide pristine 
information on the quark-gluon plasma. In addition, since the direct photon 
production in nucleon-nucleon collisions ||100|| is understood, its production 



in nucleus-nucleus collisions could provide constraints on the little-known 
gluon structure functions in the nucleus via the Compton process gq q'j. 

Data on single photons in S + Au collisions [|101|] was claimed to be under- 
stood only if the quark-gluon plasma is formed, because an overly simplified 
hadronic fireball model led to yields an order of magnitude [[102| , |103|| above 
the data. The problem with this interpretation is that the fireball model is 
inconsistent with the observed rapidity and transverse momentum distribu- 
tions. 

A Large background is another problem with the direct photon measure- 
ments. Compared to dileptons, the photon background, which mainly comes 



from VTo and r] decays, is usually much higher |90|. Direct photon signatures 



are therefore rather difficult to disentangle [|104| . 



59 



1.3.4 Disoriented Chiral Condensate 

The restoration of chiral symmetry at high temperatures is one of the most 
basic features of the QCD matter. According to lattice QCD calculations 
the critical temperature for chiral restorations coincides with that for the de- 
confinement transition. Thus the chiral phase transition in QCD is expected 



from both the theory study and lattice numerical simulations pl| , |105|| . A 
possible signal for the chiral transition is the formation of a coherent pion 
field called the disoriented chiral condensate (DCC) |p.06| , |107| , |108|| . In equi- 



librium, a DCC can not grow larger than ~ l/Tc ~ Ifm ||107|| . However 



in non-equilibrium it is possible to grow larger. A large domain of DCC 
will manifest itself by the striking fluctuation in the spectrum of neutral and 
charged pions in terms of the ratio n^o/nJ^°*"'. There also exists some cosmic 



ray data ||109|| with that ratio far from the average value 1/3. 

Most studies of this phenomena are based on the linear sigma model, with 
the effective potential as 

Veff = \{^^ + a^-v^f/A-Ha , (1.139) 



where = \/2\v'^ as the sigma mass at T = if -ff = 0, The explicit chiral 
symmetry breaking term Ha takes into account the non-zero mass of the 



quarks in the real world. In the so-called quench scenario ||108|] , the initial 



(cr, 7?) field is assumed to be symmetric under chiral transformations 

cT^cr + e-7f , (1.140) 



60 



(1.141) 

with < a >=< TTj >= 0. The evolution of the field is then described by 
zero-temperature equations of motion. Long wavelength modes of the pion 
field are amplified according to 

= (Xv^ ~ , (1.142) 

because the Mexican-hat form of the effective potential makes the a = ttj = 
point unstable. In this quenched scenario the domain size could grow as large 
as 4 — 5fm [|110|| . In the annealing scenario, {a, vf) mean fields evolve almost 



synchronously with the effective potential, thus the system just oscillates 
close to the equilibrium point. Long wavelength modes are less amplified 
and domain size is therefore smaller [ 110| . 



DCC growth depends strongly on the initial condition, and it is one main 



uncertainty in present DCC studies |110, 111|. Furthermore, most studies 



are based on the linear sigma model, and this model is not dependable at the 
starting chiral symmetric phase where the temperature is above Tc. To verify 
that some pions come from the coherent source of the DCC domain, one can 
seek secondary signals such as pion pair correlations. However, Asakawa 
probed that the pion pair correlations are very small and mainly a.t 6 = 180° 
| |112| | . Theories must be improved in order to give the size of the DCC domain 
in heavy ion collisions and/or the possibility that we would observe a strong 
pion number imbalance. Experiments to detect DCC formation could be 
done at RHIC, or as proposed in |p,13| . 
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1.3.5 Directed Transverse Flow 

Transverse flow from azimuthal anisotropy was observed recently in Au + Au 



collisions at lO.SAGeV at the AGS 114 . It was first discovered in 1984 at 



much lower energies at Bevalac |[L15|| , and has been the subject of intense 
studies since then. In general, flow can provide information on the compress- 
ibility of the dense baryonic system and on the equation of state. It can also 
assist us in studying the role of the mean fleld in nuclear collisions. 

The azimuthal distribution of energy and multiplicity would be isotropic 
if there were no collective transverse flow. Transverse flow manifests itself via 
the anisotropy in the azimuthal distribution on an event-by-event basis. To 
study the anisotropy, one can Fourier-transform the distribution /(0) ||116 



= E /^e'""^ = E v^e-''^"e'^'f> . (1.143) 

i i 

For real functions, the above equation can be written as 

/(0) = E "^^ [cos{mpn) cos(?7.0) -|- sm{mpn) sin(n0)] . (1.144) 

i 

Anisotropy gives non-zero values of f„ {n > 1). Variables Vi and ipi for the 
flrst harmonic represent a shift of the center of the isotropy circle, with tpi 
being the azimuthal angle of the reaction plane, which is the plane formed by 
the beam axis and the impact parameter vector. The variable V2 represents 
the eccentricity of the ellipse in the reaction plane. 

The above analysis is complicated by the fluctuation resulting from flnite 



multiplicity in a nuclear collision |117|. Therefore the reaction plane can not 
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be exactly determined event by event, and the values of i/jn suffer from this 
uncertainty. More importantly, the above variables for the harmonics for 
an ensemble of events are described not by a constant, but by a Gaussian 
distribution with a finite width. 

The two-dimensional Gaussian would be centered at zero if there were no 
transverse flow: 

(Pw 1 ^ -, ..X 

— ; — rr = 7; — ■ 1.145 
Vndvndipn 2na:^ 

Assuming that the fluctuation changes little with transverse flow, for 

f{<t>) = T.fA<l>-<l>i) (1-146) 

i 

we have 

< - -.<Y.fl>--.<fl> ^ (1-147) 
N 

^(7„ = —(71 . (1.148) 
TT 

Note that cr^ is independent of n, and we now denote it as cr. 

With transverse flow, the Gaussian is centered at non-zero {vn,%l)^: 

(Pw 1 (vn-vn)'^ 1 Vn+Vn'^-'^VnVnCOB(ipn-'<l'n) 

-e = -e ^ , (1.149) 



VndVndt/jn 27r(T^ 27r(7^ 



where Iq{x) is the modifled Bessel function. 
Redeflning the dimensionless variable as 



-e-^/o(^) , (1.150) 



— — , (1.151) 

(7 (7 
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the new distribution becomes 

^ = ^e-(-"+^")/o(x.x„) . (1.152) 

The above distribution -^j^ resembles a Gaussian function when x„ < 1 
and then is hard to be distinguished from the finite multiphcity fluctuation. 
However, when the transverse flow is strong enough so that af„ > v^, it has 
a minimum at Xn = 0, which is a distinct signature for the transverse flow. 



From the experiment E877 at AGS |118|, clear transverse flow signal is 



observed. In mid-central events, a clear minimum in the -^^r- distribution 
is seen. The vi coefficients extracted from transverse energy anisotropy and 
multiplicity anisotropy are quite different, and the reason is that nucleons and 
pions have different flow patterns. For protons, transverse flow grows linearly 
as a function of px-, and reaches 20% of the total transverse momentum at 
Pt = IGeV. From the comparison with different transport models such as 
RQMD [|119|| and ARC ||120|| , one can study the effect of the mean field in 



the nuclear collisions. The most striking collective fiow signature of a QGP 
formation has been predicted to be a local minimum of the fiow excitation 
function at the so-called "softest point" of the nuclear matter ||121| . This 
proposed signature is being looked for at AGS. 

1.3.6 Hanbury-Brown Twiss 

Correlation functions for identical particles probe the source size and freeze- 



out time of the dense medium |122, 123 1. The two-particle correlation func 
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tion in momentum space is defined as 



C2(P1,P2) 



P{Pl,P2) 



;i.i53) 



PiPl)PiP2) 

If tlie particles are classical without dynamical correlations, the above func- 
tion is always equal to 1. For identical quanta with bosonic statistics | |124 |, 
symmetrization introduces a dependence on two momenta and 6*2(^1,^2) = 2 
without dynamical correlations. More generally, for a source described by 
seven- dimensional phase space density function S{x,'p), the semi-classical 



expression for the correlation function is given by ||123| , 125 1 

/rfWx'5(x,i?)5(x',i?)e^''(^-^') 
^ Jd^xS{x,pi)Jd^x'S{x',p2) 

where K = {pi +p2)/2, q = {Ei — Ei.pi — P2), and / d^xSi^x, 

For an azimuthally symmetric source ||126| 



S{x, p) Qc exp 



x^ + (t - tof 



usmg 

Qo = El — E2 



pi -pi 



2m 



K-q 



2RI 



2A2 



(1.154) 
= P{p)- 

(1.155) 



= P ■ q — Poutqout + Plongqiong , (1-156) 



El + E2 Ek 

where qiong is defined to be the component of q along the beam direction, qout 
is the transverse component which is in the plane with the beam direction 
and the K vector, and qside is the other orthogonal transverse direction. The 
approximation sign is valid for |g| <^ \K\, and (3side oc K^ide = by definition. 
Thus 



C2{K,q) 



1 + exp -Ri{qi + q') - R^qi - A^o 



;i.l57) 
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~ 1 + exp [-qldeRl - llutiRl + PouA^) 

~lfongiR'l + Pfong^'^) ~ '^PoutPlongQoutQlong^'^ ■ (1-158) 

In the longitudinally co-moving system (LCMS) where Kiong = 0, it be- 



comes 



C,{K, q) = l + exp -gLe^i - lltiRl + PouA') - iLA ■ (1-159) 



Experimentally one usually uses the parameterization [ |127 



C2{K, q) — 1 + A exp —q^sideRiide ~ QlutRlut ~ QlongRlong 



;i.i6o) 



Thus one can obtain information about the space-time dimension and ex- 
pansion velocity of the emitting source from the fitted radius parameters. 
The chaoticity parameter A is different from 1 due to other correlations and 
resonance decays. For charged particle correlations such as pion and kaon 
correlations, the Gamov correction due to Coulomb interactions is important 
and must be included. Data from p + P6, S + Ph and Ph + Ph collisions ||128| 
show that the fitted radii grow with the size of the system, and are generally 
larger than the projectile size. Pions have larger radii and smaller chaotic- 
ity than kaons due to more resonance decays. Because of collective flow 
phenomena associated with transverse expansion, the above simple formulae 
( |1.158| - pri60D are modified. The fitted radii do not correspond directly to the 
source size but must be interpreted in terms of a detailed transport model. 
Pion interferometry can nevertheless serve as a powerful signal of the QGP 
formation. Given a rapid crossover region around Tc, the speed of sound 
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acquires a local minimum near Tc- As matter cools and crosses that mixed 
phase region, the minimum of causes the transverse expansion to stall for 
some time. This leads to the time delay signature of the phase transition 
that could be seen in detailed systematics of the qout, Qside dependences of the 



function C2 (gout, g«de) [pl 



1.4 Open Charm 

One of the main uncertainties in heavy ion physics is the initial conditions of 
the dense parton system. As shown in section |r.2.4] , for example, the initial 
momentum distribution, energy density and chemical composition determine 
the fate of the QGP. If the assumed initial conditions in eq. ( 1.121|) 



are true, we may not have fully equilibrated quark-gluon plasma at RHIC. 
Recently it is found that minijets lead to initial conditions characterized 
by large fluctuations of the local energy density and of the collective flow 
field. The "hot spots" ||129|| corresponding to higher local energy density 



may change the usual picture from cylindrically symmetric, homogeneous 
quark-gluon plasma to turbulent, volcanic plasma. Certain signals of QGP 
may be sensitive to this change of the initial condition. Initial conditions 
are themselves determined by the parton structure functions from the two 
colliding nuclei, mostly by the gluon structure function. Parton distributions 



in nuclei are depleted due to nuclear shadowing effects ||130| , |131|| . As we go 



to higher energies in search of the quark-gluon plasma, we can see softer 
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(small-x) partons deeper in the shadowing regime. It is therefore important 
to understand nuclear parton distributions in order to calculate reliably all 
the other observables in A + A collisions. 

Open charm production is well suited to measure the initial gluon struc- 
ture function and its nuclear modifications. At collider energies ^/s > 200 



AGeV the initial minijet plasma is mostly gluonic ^ with a quark 
content far below its chemical equilibrium value. Because open charm is 
produced mainly through gluon fusion, it can provide information on both 
the incoming gluon structure functions and the later-formed dense minijet 
plasma. 



This thesis on open charm production ||132|| in high energy nuclear col- 



lisions is motivated by two recent studies which predicted widely 

different rates. In the pre-equilibrium contribution was found to be al- 



most equal to the initial gluon fusion rate, similar to the suggestion of the 
charm enhancement from thermal productions in the hot-glue scenario |^ 



In 1^^, a more provocative claim was made that open charm may even be 
enhanced by over an order of magnitude above the initial perturbative QCD 
(pQCD) rate. 

The proposed method to measure charm meson (D/D) cross section at 
the Relativistic Heavy ion Collider (RHIC) is efi coincidence measurements, 
or dilepton measurements. Traditionally, dileptons from charm decay are 
regarded as background, while thermal dileptons or Drell-Yan dileptons are 
considered as signals. Dileptons from Drell-Yan process probe the quark 
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degrees of freedom, and thermal dileptons tell us about characteristics such 
as temperature and phase transition in the dense plasma. However, those 
signals will probably be overshadowed by the large background of dileptons 
from charm meson (D/D) decay in Au + Au collisions at RHIC and 
it is easier to measure the dileptons from open charm decay. In the QCD 
factorization limit, open charm mainly comes from initial gluon fusion, and 
it thus depends strongly on the gluon distribution in the nuclei. 

The thesis discussed at length in the next chapters consists of two parts. 
First, it is shown that the pre-equilibrium charm production is small com- 



pared to the initial fusion rate [|132|, EM. Second, it is shown that this domi 



nance of the initial gluon fusion mechanism can be used to probe the gluon 



shadowing in nuclei |133|. While the first result of the thesis work is pes- 



simistic in the sense that it indicates, contrary to an earlier claim [HI], that 



open charm is not a good probe of the QGP formation, the second result 
is optimistic as it provides a new tool to measure the gluon shadowing via 
p + A reactions. 

Perturbative QCD analysis ||131| , |134|| predicts that the gluon structure 



will be depleted at x < 0.01 due to gluon recombination processes. Shad- 
owing of the quark nuclear structure functions, 5^(2;, Q^), is well estab- 
lished from deep inelastic iA iX reactions ||130| . For heavy nuclei, the 
quark structure functions are shadowed by a factor, ^ 0.1,(5^) ~ 

1.1 — 0.1 A^^^, nearly independent of Q^. It is far less known about the gluon 
shadowing. In high energy nuclear collisions, the gluon nuclear structure is 



69 



of central importance since it controls the rate of minijet production that de- 
termines the total entropy produced at RHIC and higher energies ||135|| . The 
understanding of the gluon shadowing is therefore of fundamental interest in 
this field. 

In Chapter |^ the dependence of the direct pQCD rates for charm produc- 
tion on structure functions, momentum scale is reviewed, and we compare 
our results with a previous one |Q. Chapter ^ presents the results of part 
one of this thesis on the pre-equilibrium charm production. The important 
role of space-time and momentum correlations in suppressing that compo- 
nent of the charm production is emphasized. Chapter ^ contains the second 
part of this thesis. The dilepton spectra from open charm decay in p + Au 
collision at 200 GeV/A are calculated and its sensitivity to the nuclear shad- 
owing is quantitatively estimated. We pay particular attention to possible 
backgrounds of the dilepton spectra for the PHENIX detector geometry at 
RHIC. Discussion and summary are given at the end of each chapter. Chap- 
ter H presents an outlook of some open questions on charm productions. 
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Chapter 2 

Initial Charm Production 



In this chapter we will study how the direct pQCD rates of charm production 
depend on structure functions, momentum scale Q^, and K-factor. We will 



then compare our results with results from the Parton Cascade Model |]4^ . 



2.1 Production Mechanisms 



Heavy quark production in pp reactions was calculated long ago in ||136 
including both fusion and heavy flavor excitation processes in the leading 
order pQCD. It was proposed that the flavor excitation processes were dom- 
inant at high energies because a small exchange can easily liberate any 
charm component in the nucleon while gluon fusion was suppressed because 



> 4M^. In the Parton Cascade Model |^^, both mechanisms are incorpo 



rated to calculate s, c, b quark production in nuclear collisions. Results from 
this model suggested that the flavor excitation of the charm quark of nuclear 
structure functions would be the dominant source of charm production in 
nuclear collisions as well. 
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However, it is pointed out [|137|| that the original flavor excitation rates in 



| 136|| were too high in the x/ ~ region due to neglected interference with 



other pQCD amplitudes of the same order. When all diagrams were added 
together, a large destructive interference was found to suppress the flavor 
excitation rates by powers of A/Mg, where A ~ SOOMeV is a typical QCD 
scale, and Mg is the heavy quark mass. The suppression factor appears in 
the process g + c(c) where the charm is evolved from the structure functions 
using perturbative QCD, as also shown in [138|. 



We note that there is a possible non-perturbative charm component (in- 
trinsic charm) in the nucleon. However, there are experimental constraints 
on the amount of this non-perturbative charm component [|139| , |140|| . The 



total contribution of the intrinsic charm was shown in ||141| , |142|| to be small 



(about 10%) in the midrapidity region where most of the charm is made. 
Although the contribution of the intrinsic charm component appears impor- 
tant at large xj, its contribution to the total cross section is small and well 
within the uncertainties from other sources which we will discuss. 

2.2 Cross Section 

We therefore only include fusion processes for the parton level cross sections 
as in p2|. The cross sections are the following jl36[|: 



da = dxidx2fi{xi,Q'^)f2{x2,Q'^)d-{s) , (2.1) 

where s = X1X2S. 
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Let xl = x\ — and F = xf as the fractional-momentum parton 
distribution function, then 



da = dsdx l 
= dsdxL 



Fi{xi,Q^)F2{x2,Q^)a{s) 

XiX2{Xi + X2)s 

Fi{xi,Q^)F2{x2,Q^)ais) 



sJx\ + As/ s 



where 



(2.2) 



Xl = {xl + Jxl + 4s/s)/2, X2 = {-XL + Jxl + 4s/s)/2 



(2.3) 



By assuming that only depends on s, we can explicitly obtain the 



pQCD differential cross sections for a + 6 —> cc + X as the following [ 136 | 

87ra2(g2A 



' qq^cc 



(2.4) 



3s 



'gg^cc 



-(7 + 



31M2 1 



AMI ^ 



log 



1 + X 

1-x 



(2.5) 



where x = y 1 — 4M2 / s 



2.3 Parton Structure Functions, Momentum Scale, and 
K- Fact or 



For the production in pp collisions, we use the light quark and gluon struc- 
ture functions from Gliick et al. ||143|| and Duke-Owens |[144|| for comparison. 
In Figure |2.1| we compare the parton structure functions from these two pa- 



rameterizations together with MRSA ||145|| parameterization. The MRSA set 



parton distribution functions 
Q'=lOGeV' 
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Figure 2.1: Parton momentum distribution functions for gluons and u sea- 
quarks as a function of Bjorken x for GRV-HO, DOl and MRSA parametriza- 
tions. 
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Parton distribution functions 


A(GeV) 


GRV-LO set 


0.25 


GRV-HO set 


0.20 


Duke-Owens set l(DOl) 


0.20 


Duke-Owens set 2(D02) 


0.40 



Table 2.1: Lambda scales in different parton distribution functions 



of parton distributions is very close to the latest parameterization CTEQ4M 



We consider the following two choices for the scale in the coupling 
constant a^Q^) = 12n/ [(33 - 2n/) \og{Q^/A^)]. The scale in the parton 
distribution functions is taken to be the same. 

1. for gg — >• cc, = s/2; for qq cc, = s. {Q"^ choice-1) 

2. for both gg ^ cc ^ qq ^ cc, = s. (Q^ choice-2) 

We take Uf = A for charm quark production and Uf = 5 for bottom quark 
production. The QCD scale A depends on the choice of parton distribution 



functions and is given in table |2]l 



To incorporate approximately the next-to-leading-order(NLO) corrections 
to the above rates, we multiply the leading order results by a K-factor. In gen- 
eral, the K-factor depends on the choice of parton distribution functions, the 
center of mass energy of the collision, and the type of the projectile and target 
particles. Calculations to order 0{al) for the subprocesses were carried out 
1471 , |148|| , and afterwards the calculations to order 0{al) for p + p collisions 
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were made pi pOl. For DOl, = 1.5 GeV, = 4^2, P^^^ = 100 - 1000 



GeV, the K-factor for p + p collisions ||149|| was found to range from 2.85 to 
4.1. There are other NLO calculations of charm production cross sections 
| 151| , |152|| , which reveal the dependence of the NLO K-factor on the trans- 
verse momentum and rapidity of the initially produced charm. 



2.4 Comparison with Low-Energy Experimental pp Data 



In Figure |2.2| we compare the so-calculated charm cross section to the limited 
data on inclusive cc production in pp collisions. 



The NA34 data for cTcharm is taken directly from [153]. The values for the 
other data lines are computed from D-meson cross sections according to the 
argument in |p.54|| by using the published experiment results ||155| , |156| , |157| , 
|158|| . By assuming that Ds/Dg production accounts for 8% of the total cross 
section, hidden charm accounts for 1%, and AcD/AcD production is smaller 
than 6.1 fib at NA27, the conversion is done as the following: 



a{pp ccX) 



a{D/D) 



1 + 



a(< Q.l] 



(2.6) 



2(0.92 ±0.04)0.99 V" ' 30.2 ±2.2^ 

We then take a = 6.1(0.) as the upper (lower) limit. Including the already- 
converted NA34/2 data, we have the following cross sections for cc in ta- 
ble ^ 

Only statistical error is included in the above calculations. Earlier exper- 
iment results |159| show larger uncertainties among different experiments. 



In Figure 2.2, the solid line is our result with = 1.3 GeV, K = 3, Q 



76 



"1 r 



"1 1 1 r 



NA34 



NA27 







. E743 



E789 



/ / 



// 



/' 



m^=1.3GeV, GRV-HO, Q -choice 1, K=3 
m^=1.3GeV, GRV-HO, Q'-choice 2, K=4 
m^=1 .5GeV, DO set 1 , Q'=4m/, K=3 
m =1.4GeV, GRV-HO, O'-choice 1 , K=3 



100 200 300 400 500 600 700 800 900 1000 

Piab(GeV) 



Figure 2.2: Cross sections for pp 



ccX as a function of Piab- 
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Experiment 




published original data(/x6) 


converted accit-ib) 


NA27 


400 


^ (j{D/D) = 30.2 ±2.2 


14.7-22.1 


E743 


800 


^ a{D/D) = 48l^° 


21.0-40.5 


NA34/2 


200 




14.1 ±9.3 


E653 


800 


a{D/D) = 76 ± 10 


34.7- 60.1 


E789 


800 


^ cj{D/D) ~ 35.4 ± 1.3 


17.9-25.7 



Table 2.2: cc cross sections from low energy charmed meson data 



choice- 1 and GRV-HO set. As a consistency check, we also plot the long- 
dashed curve using the same parameters as in Figure of [ |149|| (i.e. DOl, 
Mc = 1.5GeV, = 4M^) using constant 7^ = 3 for simplicity. Both the 
solid curve and the dashed curve fit the low energy data reasonably well, so 
we use these two parametrizations for the following high energy calculations 
in this chapter. Comparison of the solid and dot-dashed curves shows the 
strong dependence on the assumed charm quark mass for the GRV-HO set. 
Comparing the solid and dashed curve we see that different choices for the 
scale can be compensated for by shifts in the K factor. These results, together 
with the large uncertainty of data, emphasize the need to measure pp and 
pA reactions in order to control uncertainties in the initial charm production 
rate, so that charm productions in AA can be properly calculated. 

As a further check of the parameters we compare charmed hadron Xf 
results in Figure with 400 GeV pp data |[L55|| . The solid curve represents 



the result da/dxj using the first parameterization. The dashed curve repre- 
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sents the result using the second parameterization. Both curves assume a S 
function charm fragmentation function. The reahstic fragmentation function 



used in [|141|| lowers the curves slightly and reveals the true high-x/ intrin- 
sic charm component. In Figure |2.4| we compare bb production. Here we 
take Mh = 4.75 GeV as in ||15CI|| , with K = 3,nf = 5. The data point at 
^/S = 630GeV is from [|T6g]: a{pp ^ b + X) = 19.3 ± 7{exp.) ± 9{th.)fxb, 
and only the experimental error is indicated in Figure p.4| . The dashed cross 
at = 1.8 Tev is obtained indirectly from [|15(J|| , and the error bar is only 
illustrative. At = l.STeV, our value is Al.Sfib x K = 125/i6. This is 



significantly larger than found in |150 



2.5 The A-scaling of Charm Productions 

The dependence of cross sections of charm productions on atomic number A 
in p + A collisions is usually parameterized as = a^A"'. Earlier exper- 
iments gave a ^ 0.76 ||161|| . However, those experiments did not explicitly 



reconstruct D mesons and had lower statistics. Later experiments ||162| , |158 
reconstructed D mesons and found a = 1.02 ± 0.03 ± 0.02 [|158|| , which is 
consistent with a = 1. from the Glauber geometry. We note that the value 
a = 0.76 is taken from low energy experiments, where energy conservation 
suppresses the contribution from multiple collisions. At high energies, QCD 
factorization implies that a = 1 for p + A scaling is the appropriate scal- 
ing modulo small nuclear shadowing effects. The nuclear shadowing effect 
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Figure 2.3: The production of charmed hadrons as a function of for pp 
coUisions at Piab = 400 GeV. 
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Figure 2.4: Cross sections for pp ^ bb + X vs ^/^. 
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on the A-scaling was estimated by Qiu ||163|| , and he found that the nuclear 
shadowing could only lead to a small reduction in a (a > 0.96) in low energy 
experiments ||161| ]. In the sections for initial charm productions, we neglect 
the nuclear shadowing since it is a minor effect compared with effects in our 
focus. 

2.6 Comparison with Results from Parton Cascade Model 

Now we compare our results for the rapidity density of produced cc pairs 
at Y = with results of |0. Figure |2.5| shows the energy dependence in 



the range between RHIC and LHC ( = 200 - 6300 AGeV) for Au + Au 
collisions. 

The scaling from pp results to AA is 

(dN \ "^"^ / \ 



where crfneiastic is taken from ||164|| , which is about 40.8m6 at RHIC energy and 
73.5mb at LHC energy. Glauber geometry for central high A + A collisions 
gives a = 1. 



In Figure 2.5 the solid curve is the result we obtained using the same 



parameters as for the solid curve in Figure |2.2| . The parameterization for 
the dashed curve in Figure |2.2| gives a curve higher than the solid curve by 
15% to 30%. The four long-dashed curves, curvel to curve4, are all from 
PCM calculations The top curve4 is the parton cascade model result 
for the so-called QGP formation case, including both the fusion and the 
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flavor excitation processes. That curve is much higher than our sohd curve 
because it includes the contribution from flavor excitation processes. CurveS, 
the curve with solid squares, shows the contribution to curve4 from fusion 
processes only (processes (1) and (2) in the notation of [@), and curve3 
is very close to our results. The bottom curvel is the parton model result 
extrapolated to AA from pp using the A^'"^^ scaling measured at much lower 
energies. It is lower than our solid curve by a factor of 12 at the RHIC energy 
to 5 at the LHC energy. The main source of this difference is from the A- 



dependence oip+A cross sections. The parton cascade model [44| used an A" 



scaling with a = 0.76 ||161|| instead of the value a = 1 we use from Glauber 
geometry. To demonstrate the effect in different A" scalings, we plot curve2, 
which is the parton model result scaled from pp to A A by A^^^, and it is close 
to our results. As shown by the two arrows, curve4 becomes curveS when 
the coherent cancellation of flavor excitation processes is considered, and 
curvel becomes curve2 when the high energy scaling is used. Therefore, the 



factor ~ 50 enhancement of charm production suggested in comparing 
curvel with curve4 for charm production at RHIC is a consequence of the 
inclusion of incoherent flavor excitation processes and the extrapolation from 
pp to AA via low energy scaling. The net dynamical enhancement in the 
PCM should be obtained by comparing curve2 with curveS. In that case 
Figure |2]5| leads to the expectation that the pre-equilibrium charm production 
should be comparable to the initial fusion rate. This removes the bulk of the 
discrepancy between [0] and [|^. 
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2.7 Summary 

We calculated initial charm production in nuclear collisions. We demon- 
strated the sensitive dependence on the choice of structure functions, the 
scale, and the K-factor. The parameters were fixed by fitting the lim- 
ited available experimental data at lower energies. We emphasized the need 
for new measurements of pp and pA charm production to reduce the present 
large theoretical uncertainties. We argued that the copious charm production 
predicted in was mainly due to the neglect of the coherent suppression 



of flavor excitation processes and the low energy pp to AA scaling. Our cal- 
culated initial charm yields are close to those computed in |^ and to the 
curve 2 in Figure |2.5| from . 
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Figure 2.5: Comparison with Parton Cascade Model results on 
{dNcc/ dY)Y=Q- the rapidity density of charm and anti-charm pairs for 
Au + Au collisions as a function of y/S/A. 
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Chapter 3 

Pre-equilibrium Charm 
Production 



In this chapter we consider the pre-equihbrium contribution to the charm 
yield in A + A colhsions. It is produced through final state interactions 
between partons in the dense minijet plasma. Here we only calculate the 
dominant contribution from minijet gluon fusion. 



3.1 Spectrum of Minijets 



The spectrum of minijet gluons in leading-order follows from [|165 



da 


9 7ral 


dt gg^gg 


2 P 


da 




dt gq^gq 


s2 



^ ut us st 

s2 p ijp' 



Au^ + s^ + s^ 
+ 



(3.1) 



9 us t2 

The term minijets refers to unresolved jets at a scale j9j_ > Pxcut- The 
inclusive cross section to produce minijets is given by 



da 

dydp\ 



j dy,x^^X2f2^il + 2^3 + A) 



(3.2) 
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where /i is the incident parton distribution evaluated at xi = p_L(e^+e^^)/ y/s 
at a scale = p±^. The light-cone coordinates of the initial and final partons 
are: 



Pi 

P2 
P3 



2xiPo,0,0 
0, 2x2Po, 



m^e^^, m±e —pi 



(3.3) 



The observed parton has momentum 



P 



(3.4) 



The subprocess Mandelstam variables are s, t, and u. For the calculation of 
minijet gluon fusion process in the following Chapter ^.2|, we choose = 



s. As in [£2], we use DOl as the proton structure functions and K = 2, 
Mc = l.SGeV for the minijet production. The nuclear shadowing effect is 



taken from ||166|| . In Figure we show the resulting transverse momentum 
distribution of mid-rapidity minijet gluons at ^/s = 200 AGeV in open circles. 
We call this distribution the hard distribution since it only includes the hard 
gluons with p±_ > p±cut = 2GeV . The solid line in Figure O] is the output of 



the Monte Carlo calculation via the HIJING model that includes initial 
and final state radiation, and therefore fills the soft gluon region. The dashed 
line is our fit 0.06e~^'^^^^, which is discussed below. In Figure |3.2| we show 
the rapidity distribution of minijet gluons. The circles represent our result 



from the initial production. The dashed line represents the fit cos 
where dN/ dy scale is arbitrary. 
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Figure 3.1: The minijet gluon distribution A '^1'^ [dN / dydp]^ 
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Figure 3.2: The minijet gluon rapidity distribution dN/dy. 
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For convenience we have parameterized the minijet spectrum (both soft 
and hard gluons) from the Monte Carlo calculation as the following: 
dN 



, ,^ - g(p±)p(y)A^/^ = 0.06e-^-25f^ cos 
dydp± 

with \y\ < 3.7 . 



A""/' , (3.5) 



In the following, we call this parameterized distribution the soft+hard distri- 
bution. The soft+hard, hard, and Monte Carlo distributions are very close 
to each other in the semi-hard p± > 2 GeV region at y = 0, as seen in Fig- 



ure |3.1| . The parameterized distribution falls underneath the Monte Carlo 
result in the region p_L < 1 GeV. However, the soft component is strongly 
model-dependent as it requires the furthest extrapolation from the pQCD 
hard domain. The Hijing yield in that region is due to initial and final state 
radiation. Other contributions in this soft domain from coherent strings are 



possible [^. While most of the following results are obtained with the sim- 
ple parameterization above, we will check the sensitivity to variations of the 
soft component as well. We also note that at bigger rapidity the p± spectrum 
falls more rapidly. The above parameterization does not include that prop- 
erty. However, this property only lowers the high p± tail, and only slightly 
changes the low p± part and the total number of the pre-equilibrium charm. 



3.2 Two Extreme t] - y Correlations 

In ideal Bjorken dynamics, the space-time rapidity 

7^ = l/2log[{t + z)/{t-z)] , 



(3.6) 
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and the true momentum rapidity 

y = l/2\og[{E + p,)/{E-p,)] (3.7) 

are assumed to be perfectly correlated. This is referred to as the inside- 
outside picture, and the phase-space distribution function in Bjorken corre- 
lation case has the form 

t[X,p,t)Bj — ^ 



T'7iRj^p± dydp± 
X 6{r]-y)Q{T-T,)Q{Tf-T) , (3.8) 

where = 0.1/m/c is the minijet formation time. Tf ^ 1.7 fm/c as the 
proper time when the energy density of the pre-equilibrium minijets falls by 
an order of magnitude to ~ 2GeV / fm? due to rapid longitudinal expansion, 
and that is when we terminate the pre-equilibrium stage. 
The phase space distribution is normalized such that 

r F{x,p,t)Bjd^x ^d^N ^ 1 dN 
J (27r)3 ~ d^p ~ E dydp-1 ' ^ ' ' 

In this chapter we study the pre-equilibrium charm production at = p2 



'd^N\ f A f 1 d^pid^p2d^p' 

E—— = d X 



d^pjy^, J J 32(27r)8 uj,uj2E' 

X F{x,pi,t) F{x,p2,t) |M|2 5(4) (^pA«) . (3.10) 



Denoting 



dN/dydp~l = g{y,p±) , 

pli = p±i(cos0i,sin0i,O) , (3.11) 
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the ideal rj — y correlation leads to 



E 



g{v,P±i)9iv,P±2KJ:E)\M\^ 



X J dridp_Lidp±2d(f)id(f)2- 



E' 

ln(r//r,) 



32{27if7rRl 



dridp±2d(pid(f)2 



^ 9{v,P±i,o)9{v,P±2)\M\'^ 

P±2 [1 — COS(01 — 02)] — {E cosh Tj — p COS (pi) 

In deriving the above, we have used kinematic relations 

E' = {p^^+p±2)coshr]- E , (3.13) 

s(j:e) 5(p±i-p±i,o) 



E' p±2 [1 — cos(0i — 02)] — {E coshr/ — p cos0i) 



(3.14) 



^ p±2(-^cosh?7 -pCOS02) ,g 

pj_2 [1 — cos(0i — 02)] — (E coshr] — p COS0i) 



In Figure we show the numerical results of the above integral in equa- 
tion (^J^). 

The solid line is the p^-distribution for the initial charm production 
from Chapter |^. The curve labeled with solid diamonds represents the pre- 
equilibrium contribution including both the soft {p± < 2 GeV) and hard 
{p± > 2 GeV) components of the minijet gluons. The curve labeled with 
hollow diamonds represents the pre-equilibrium contribution including only 
the hard component. We see that the pre-equilibrium contribution in this 
strongly correlated case is totally negligible. This result is similar to the 
thermal charm production contribution calculated in ^2\, except that, in 
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Pj (GeV) 



Figure 3.3: The distribution [Ed^N / d^p)y^Q of charm quark production as 
a function of p±_ using 5{ri — |/) -correlation . 
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our case, the curve extends to higher because of the broader initial mini- 
jet distribution in p±. 

Another extreme case, opposite to the ideal Bjorken picture, was con- 
sidered in ||6^. In that uncorrelated r] — y case, the gluon distribution was 



assumed to be completely uncorrelated as in an ideal thermal fireball. This 
assumption leads to 

(3.16) 

p V ayap± 

If one assumes a fixed volume V = TinR^, then J dt Tf — Ti, and 
Then from equation (|3.1CI| ), we have 



' d'N\ _ Up) fd'x _ Tfln 



where 



j^^-^ ^ 1 dyidy2dp^2d<Pid<p2 9iyi,P±ifl)9iy2,P±2)\M\'^ 
J cosh ?/i cosh ^2 D{yi,y2,p±2,4>u4>2,p) 

E' D{yi,y2,p±2,(pu(l)2,p) 

p±^{E coshy2-pcos(f)2) 
P^ho = I T ' 3.21 

Diyuy2,P±2, </'l,</>2,P) = P±2 [C0sh(?/i - 7/2) - COs(</)i - (j)2)] 

— {E coshyi — pcoscpi) . (3.22) 



In Figure ^.4| we show the pre-equilibrium charm production for this uncor- 
related case. The solid curve represents the initial charm production. The 
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curve labeled with solid circles represents the pre-equilibrium contribution 
including both the soft {p± < 2 GeV) and hard {pj_ > 2 GeV) components. 
It is much larger than the Bjorken-correlation case, and is comparable with 



the initial charm yield, as shown in Figure |3.4| . This is similar to the result 



m 



62| in which the pre-equilibrium charm production has almost the same 
magnitude and p^-shape as the initial charm. The curve labeled with hollow 
circles represents the pre-equilibrium contribution including only the hard 
component. 

3.3 Minimally-correlated r] — y 

From the above results we conclude that the pre-equilibrium charm produc- 
tion is very sensitive to the f] — y correlation of the minijets. We consider 
here the simplest source of 1] — y correlations resulting from the minimal 
geometrical spread in initial production points required by the uncertainty 
principle. This type of correlations are included in the parton cascade model 



and discussed in [ |167|| . The phase space distribution function including such 



minimal correlation has the form 



dN d{Trr 



At 

3. 



X pQ(xQ,tQ)6{x — xq — vAt)d xodtQ . (3.23) 

The integration is over the space-time coordinates {xo,to) of the production 
points of the gluons. These points are distributed according to a normalized 
density po{xo,tQ). The delta function arises to take into account the free 
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Figure 3.4: The distribution {E(PN/cPp)y^Q of charm quark production as 
a function of p± for the non-correlation case. 
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streaming of the partons from the production point, with velocity v = p/E, 
where E = p±cosh.y and Pz = p±smh.y. The theta function defines what 
we mean by pre-equihbrium. The proper time when the pre-equihbrium 
fusion is terminated is r^ax, which is determined below in Figure |3.6| . The 
theta function insures that only those gluons with proper time less than Tmax 
contribute to pre-equilibrium charm productions. 

The formation physics is included via the Lorentzian formation factor 
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[l + {tf{p)/At)r' , (3.24) 
where At = t — to is the elapsed time, and the formation time is given by 

tf(p) ^ coshy ^'^^^^ f/m) , (3.25) 

P± 

We note that, at the early time, the above formation factor more accurately 
describes the interference phenomena suppressing production than the con- 
ventionally assumed factor 

e[At-tf{p)] , (3.26) 

We assume 

Jpoixo,to)(fxodto = l ■ (3.27) 
In this case the normalization factor is 

AT = ^ , (3.28) 
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so that 

lim — — / F{x,p,t)Mind^x = —- . 3.29 



As discussed in ||167|| , the production points are spread along the beam axis 
according to the uncertainty principle by an amount 6z = d ~ /i/p_L, since 
the dominant parton interaction leading to a ?/ = parton with final p± has 
an initial longitudinal momentum xPq ~ p^. We take as a particular model 

d=—(fm) . (3.30) 

Clearly this is only a rough guess, but it allows us at least to investigate the 
sensitivity of the results to a particular rj — y correlation that results from this 
spatial spreading of the production points. We emphasize that it is precisely 
the uncertainty of the initial space-time formation physics that leads us to 
investigate the possibility of open charm production as an experimental probe 
of this physics. 

Given the above assumption we take 

Po{xo,to) = —-^6{to) — j=r— , (3.31) 



where d is the mean spread for gluons depending on p±_ from above. This 
distribution only spreads out the production points along the beam axis. A 
more realistic treatment would also include a spread in the time coordinate. 
Neglecting transverse expansion, we obtain finally 



pdy dpi 1 + (^^^)2 ■ ^ ' ' 
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Denote 



ai = tanh yi 
a2 = tanh 1/2 



0.2 
^0.2 



/2 



(3.33) 



then after integration over z, we have the final expression as the following: 



E- 



PX1,0 P±2| /\f|2 
^^1,0 0.2 0.2 l^^-' I 



4p2 _/ dy2dyidp^2d<p2d(l)i 

g{yuP±ifl)9{y2,P±2) 



P±2 COsh|/i COsh?/2 (-E COsh?/2 — P COS02) 



V57T&2 



X / dt 

10 



g i/f-i+i/ba 



T _|_ / 0.2coshyi \2 



1 _l_ / 0.2coshy2 \2 



(3.34) 



In the above 



tf = Tmax min (cosh yi, cosh 7/2) 



(3.35) 



and p±ifi is the same as in equation ( ^.22] ). Note that by using the unit 
GeV for momentum and unit fm for time, the expression m 



equations ( |3.12| ), ( |3.18| ), and ( p.34| ) has the dimension GeV 



m 



We 



require a factor {hc}'^ ~ (0.26*6^/7?^)^ to convert it to the dimension GeV ^, 



which we have used in Figure p73| , Figure p^, Figure pT5| and Figure |3T7 . 

The numerical results of the above equation 3.34 are shown in Figure 
The curve labeled with solid squares includes both components while that 
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labeled with hollow squares includes only the fusion of hard gluons. The curve 
labeled with diamonds represents the contribution from fusion of soft gluons 
both with pj_ < 2GeV. This shows that the pre-equilibrium contribution 
comes mainly from the fusion of soft and hard gluons. 

In Figure |3.6| the energy density at z = as a function of proper time is 
shown assuming minimal correlation and Lorentzian formation probability. 
The solid curve includes both soft and hard components while the dashed 
curve is calculated using the hard distribution and includes only the hard 
component. The energy density increases and reaches maximum at the time 
about 0.1/m/c, then the energy density decreases linearly to ~ 2GeV / fm? 
at ~ 1.7/m/c for soft+hard distribution. We choose the above time as the 



3.4 Why Is the Pre-equiUbrium Charm Yield So Small?: 
An Estimate 

To understand why the pre-equilibrium charm yield is so small compared 
to the initial yield as found through tedious numerical calculations in the 
previous chapter, we consider here the calculation of the total number of 
pre-equilibrium charm pairs. The expression for that number is given by 

^ = / / ^^^nSJumS.Msa{s) , (3.36) 

where a{s) is the integrated cross section for the process gg cc, see equa- 
tion (p.5|) . Our main strategy is to estimate the mean difference between the 
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Pj (GeV) 



Figure 3.5: The distribution [Ed^N / d^p)y^Q of charm quark production as 
a function of using minimal r] — y correlation. 
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Figure 3.6: The energy density at z = as a function of proper time. 
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two gluon rapidities, then, from the kinematical constraint on charm pro- 
duction (s > 4:M^), estimate the effective lower cutoff for p± of the minijet 
gluons. Thus we separate the p± integrals from the rapidity integrals, and 
have a rough estimate for the total number of charm pairs. 
For the fireball case, 

Fix,P,tU = -^y^^ ' (3-37) 
s = 2p_LiP±2 [cosh(?/i - 1/2) - cos(0i - 02)] • (3.38) 
For the Bjorken case, 

F(l,Ms, = <>lJ^'±glo(r - rmrf - r) . (3.39) 

Pj_ dydp± TTiRj^ 

s = 2pxiP±2 [1 - cos(0i - 02)] . (3.40) 

For all the cases, we use the fit to the gluon distribution given by equa- 
tion (|3.tj|), where 



9{Pl) = ae-^P^ = OMe-^-^^^^ . (3.41) 

Therefore, 

ArcR^ J coshyi J cosh?/2 

X J dp±igip±i) J dp±2giP±2) J #1 J d(j)2sa{s) , 

dp±igip±i) J dp±2giP±2) J #1 J d(j)2sa{s) . (3.42) 



Nbj = I dv[piv)] 



X 
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The dominant contribution is coming from the vicinity of the production 
threshold where s = 4M^ = 9GeV^ [|136|| . From this we make the following 



rough estimates: 

sa{s) ~ q;^(s) ~ 0.06 , 

dp±g{p±) ~ dp^g{p^) ~ , (3.43) 

where pc is the effective cutoff value for p±i and p±2 from the requirement 
s > 4M2. 

For the fireball case, 



(C0sh(|/i - 7/2)) = -coshyicosny, Q ^ ^3^^^^ 

J t>± ^■'cosnyi cosnj/2 

Since the minijet p± spectrum is dropping almost exponentially, the produc- 
tion heavily favors the smaller cutoff pc, so the mean value of cos(0i — ^2) is 
most likely to be negative. We take 

(cos(0i - 02)) 0.5 . (3.45) 

Then 

sr^9pl , (3.46) 
so the effective cutoff for the fireball case is 

Pc ~ l.OGeV . (3.47) 
On the other hand, for the Bjorken case, 

yi=y2=V , ^ s ~ , 

=^ ~ 1.73GeV . (3.48) 
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Using the same values as in Chapter |3.2| : Xj = 0.1/m, Tf = 1.0/m for fireball 
case, rj. = 1.7 fm for Bjorken case, and 

dyi — ~ 2.8 , 

cosh yi 

J dri [p{ri)f ~ 4.9 . (3.49) 
We then have the estimate for the total number of the pre-equilibrium charm: 
N,, ~ l/i^^ 2.8^(^)^(2vr)^[ga(g)]^3.5 , 



Nb, ~ -^^2 4.9(^^)2(27r)^ [sa{s)] ~ 0.098 . (3.50) 



Therefore we estimate Npb/NBj ~ 35, in rough agreement with the detailed 
numerics. We see that the large increase going from the Bjorken case to 
the fireball case results mainly from the different p± cutoff. In the uncor- 
related fireball case, one allows particles with different rapidities to interact 
with each other (see equation ( 3.44| )), thus more low pj_ gluons can take part 



in the interaction. Since the minijet p± spectrum drops almost exponen- 
tially, the fireball case produces much more pre-equilibrium charm than the 
Bjorken case (a factor of 6 increase from the smaller p± cutoff). Although 
the questionable linear proper time dependence in the fireball case also gives 
an considerable increase (about a factor of 3.5), it is not as important as the 
correlation effect. 

For the minimal correlation case, the estimate is unfortunately not as 
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straightforward. The phase space distribution function is 



F{x, p, t) 



Min 



27171 R\ he cosh y 

dN t hccoshy /'q \ 

ay ap± coshj/ p± 



and s is the same as in equation ( |3.38| ). In the above distribution function 
we choose to use the ^-function for the formation-time effect. 

Using equation ( p. 361) and after the integration over z, we have 



Nuin = ^^^2 ^ J dP±i9{P±i) J dp±29{p±2) jd<Pi 

X f,y,AfL f,y,l^Sa(s) r~M^^^ .(3.52) 

J coshyiJ cosh 7/2 -^W Jl/61 + 1/62 



where 



_ /^ coshyi cosh 7/2 ^ 

tmin — ftCmaX , , 

imax = i"/min(cosh?/i, cosh 2/2) , (3.53) 



and ai, a2, 61, 62 are defined the same as in equation (p.32|) 



We estimate that for the dominant part of the integral 



bi 62 \ Pc J 



he cosh y 

Pc 

tmax ~ Tf coshy , (3.54) 
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where y = {\yi\ + \y2\)/2. Now let m = t pc/ {he cosh y), then the last 3- 
dimensional integral in equation ( p.52|) without the factor sa{s) is 



cosh yi J cosh y2 2 

rfPc 

-R^ r ainh(i;i-i;2) cosh y , 2^2 

X / due 2coshyi coshj/2 ' . (3.55) 



The M-integral gives a dependence on tj which is similar to the logarith- 
mic dependence in the Bjorken case, and the exponential form in the inte- 
grand forces the spread yi — y2 to be small. Numerically, by taking r/ ~ 
1.7/m/c , pc ~ 2.0 GeV (as the first-step value) in the w-integral, the 
above 3-dimensional integral is J ~ 19.1. When the integrand is weighed 
by cosh(yi — 1/2), the integral is ~ 23.3. So 

(cosh(yi - 2/2)) ~ 23.3/19.1 ~ 1.22 , 

^ s ~ 3.44p2 , ^ ~ 1.62GeV . (3.56) 

Note that the above determined value of Pc is insensitive to the first-step Pc 

value we tried in the w-integral. 

Therefore for the total pre-equilibrium charm number, 

(ficj'^A^^'^ r r f f 

Nmiu ~ 47ri^2 ^ J dp±i9{p±i) J dp±29{P±2) J d<f)i J d(j)2 J sa{s) 

Therefore NMin/Nsj ~ 1. From the above estimate we can see that although 
the minimally-correlated case allows particles with different rapidities to in- 
teract, the dominant contribution still comes from the region where the two 
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gluons have almost the same rapidity, thus there is no sizeable enhancement 
in the pre-equilibrium charm yield. The Minimally- Correlated case is very 
much like the Bjorken case in that the dominant contribution comes from 
Ui ~ y2 region. 

As a comparison to the above rough estimates in this chapter, the nu- 
merical integration gives Npi^ = 3.8, Nbj = 0.093, and NMin = 0.078, so 

NMin/NBj ~ 80%. 



3.5 Discussion 

In order to see the formation-time effect, in Figure p.7| we compare the results 
from two different formation-time assumptions. The solid curve is obtained 
using the Lorentzian form in equation (|3.24|) , and the dashed curve using the 



theta function form in equation ( |3.26|) . The result from this 6'-function is 



about 10% higher at p± = OGeV, and 10% lower at p± = 9GeV. The lack of 
sensitivity to the formation-time physics is due to the relative large p± for 
the gluon minijets in the charm production process. There would be more 
sensitivity had the production been dominated by low p± components. 

We also see that for the soft+hard distribution, soft gluons significantly 
increase the pre-equilibrium charm production in both low-p^ and high-px 
region, with the largest increase in low-p^ region. It is interesting to identify 
where the enhancement comes from. In Figure ^.5|, the curve with diamonds 



shows the contribution from the fusion of soft gluons both with p± < 2 GeV, 
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Pj (GeV) 



Figure 3.7: The distribution [Ed^N / d^p)y^Q of charm quark production 
using different formation-time probabihty distributions. 
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and the curve with hollow squares shows the contribution from the fusion 
of hard gluons both with pj_ > 2 GeV. These two curves are both very 
low compared with the curve calculated from the soft+hard distribution. 
Therefore the enhancement going from the hard distribution to the soft+hard 
distribution is mainly due to the fusion of hard and soft minijet gluons. 

We have noted before that our fit for the minijet gluon spectrum falls 
below the Monte Carlo result from HIJING calculation. We can fit the soft 
gluons from HIJING better by using 0.265e-^-^P^ for p± e (0, 1.1) GeV, 
and use the old fit 0.066"^-^^^-^ for higher-p^ gluons. This new fit gives us 
more very soft gluons. Using the new fit, we have done the calculation for 
minimally-correlated case, and the result is different only by less than 10%. 
We therefore conclude that super-soft gluons are not very important for the 
pre-equilibrium charm production. 

3.6 Summary 

For the pre-equilibrium charm production, we studied the effect of correla- 
tions between the rapidity y and space-time rapidity rj of minijet gluons. In 
the ideal Bjorken-correlated case, where rj — yi — 1/2, the pre-equilibrium 
charm production is negligible compared with the yield due to initial gluon 
fusion. For the opposite extreme fireball case, corresponding to uncorrelated 
y and r), the pre-equilibrium charm production is almost a factor of 50 larger 
than in the Bjorken-correlated case and is comparable with the initial charm 
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yield |62|. With the estimates of the total pre-equilibrium charm number, 
we found that the difference comes mainly from the rj — y correlation. There- 
fore, the pre-equilibrium charm production is very sensitive to the [t] — y) 
correlations in the initial state. 

In order to investigate the effect of more realistic correlations that may 
exist in the initial minijet plasma, we introduced a minimal correlation model 



taking into account the uncertainty principle along the lines of ||167|| . This 
minimal correlation is similar to the ideal Bjorken correlation case and pro- 
duces negligible pre-equilibrium charm compared with the initial charm yield. 



Ill 



Chapter 4 

Gluon Shadowing In Lepton 
Pair Spectra Prom Open 
Charm Decay 



Dileptons from charm decay are probably dominant relative to the conven- 



tional signals from the plasma at RHIC such as thermal dileptons |£8[ . Here 



we offer a different point of view on the dilepton measurement |[L33|| . Instead 
of considering it merely as a background, we propose to use opposite-sign 
lepton pairs (ee,e/i and /i/i) from open charm decay in p + A collisions as a 
measure of nuclear shadowing effects. We focus here onp + A collision rather 
than A + Ato minimize the combinatorial tt, K decay backgrounds and other 
final state interaction effects. 



4.1 Introduction 

The spectrum of dileptons produced in heavy ion collisions at high energies 
has been proposed in the past to provide information about the dynami- 
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cal evolution of quark-gluon plasmas (see, e.g., the review by Ruuskannen 
169|| ). Unlike hadronic probes, this spectrum is sensitive to the earliest 



m 



moments in the evolution, when the energy densities are an order of mag- 
nitude above the QCD confinement scale (1 GeV/fm^). However, the small 
production cross sections for pairs with invariant mass above 1 GeV together 
with the large combinatorial background from decaying hadrons necessitate 
elaborate procedures to uncover the signal from the noise. The PHENIX 
detector [^, now under construction at the Relativistic Heavy ion Collider 
(RHIC), is designed to measure ee, e/x, and /i/i pairs to carry out this task. 
One of the important background sources in the few GeV mass range arises 
from semileptonic decay of charmed hadrons {DD). As shown recently in 



||88|| , the expected thermal signals in that mass range may only reach 10% 
of pairs from open charm decay. It is much more difficult to precisely mea- 
sure these signals than to measure dileptons from charm decay, and special 
kinematical cuts and precise e/i measurements are necessary to uncover the 



thermal signals p9 



There have been several attempts |44, 62 1 to take advantage of the large 



open charm background as a probe of the evolving gluon density. Most 



mid-rapidity charm pairs are produced via gluon fusion ||142|| . Therefore, 
the dileptons from charm decay carry information about the distribution of 
primordial gluons before hadronization. In fact, the inside-outside cascade 



nature of such reactions greatly suppresses ||132|, all sources of charm 



production except the initial perturbative QCD source. Therefore, the open 
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charm background is dominated by the initial gluon fusion {gg — >■ cc) rate. 

The initial cc rate depends strongly on the nuclear gluon structure func- 
tion ||170|| , gA{x,Q^). The quantity of fundamental interest ||131| , |134|| is the 
gluon shadowing function 

D /' /n2\ 9a{x, Q"^) , . 

R./aM) = Ag^^.^Q.^ ■ (4-1) 

The point of this study is to demonstrate that the A dependence of continuum 
dilepton pairs in the few GeV mass region provides a novel probe of that 
unknown gluon structure. We also show that the required measurements of 
ee,efi, and /i/i pair yields in p + A i'^i~X at ^/s = 200 AGeV are not only 
experimentally feasible at RHIC but also that the open charm signal can be 
easily extracted via the proposed PHENIX detector. 

A key advantage of the continuum dilepton pairs from open charm decay 
over those from J/ ip decay is that it is possible to test the applicability of the 
underlying QCD dynamics at a given fixed ^/s by checking for a particular 
scaling property discussed below. In quarkonium production the mass is 
fixed and the required scaling can be checked only by varying the beam 
energy. \n p + A ^ J /ip production the required scaling was unfortunately 
found to be violated in the energy range 20 < ^/s < 40 AGeV ||171|1 , thus 
precluding a determination of Ra- Possible explanations for the breakdown 
of QCD scaling for J/ip production and its anomalous negative xp behavior 
|170| , |171|| include interaction of next-to-leading order quarkonium Fock state 
with nuclear matter nuclear and co-mover J/ip dissociation [ |172|| , and 



parton energy loss mechanisms [ 173 |. It remains an open question whether 
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the required scaling will set in at RHIC energies 60 < ^/s < 200 AGeV. 

We propose to use ee,e/i and /i/i pairs from open charm decays to probe 
the shadowing effects in nuclei. These three measurements cover three dif- 
ferent rapidity regions as a result of the PHENIX detector geometry. In 
Chapter ^ we calculate the lepton pair spectrum from open charm decay in 
p + Au at 200 GeV/A using two commonly used shadowing scenarios: the 
parameterization which does not depend on scale or parton component ||70 



and the scale-dependent, component-dependent shadowing ||166|| . The differ- 
ence between these two shadowing scenarios is significant. The ratio curve of 
dilepton dN / dMdy spectra from p + Ato those from pp has a similar shape 
to the shadowing curve. We then show that there is an approximate scaling 
among the ratio curves at different dilepton pair-mass M, which enables us 



to bring them together to the shadowing function curve. In Chapter [4.7| we 



first estimate the original backgrounds of the lepton pair spectra. The large 
number of electrons from Dalitz decay and the muons from random decay of 
pions or kaons results in large backgrounds for e/i and dilepton experiments. 
Without suppression and kinematical cuts, it is much higher than the signal 



from charm decay. In chapter |4.8| , we consider the PHENIX detector geome- 
try, a specific kinematical cut and the suppression factors of the background 
leptons. We calculate the lepton pair signal and backgrounds entering the 
detector, and the signal-to-background ratio is very promising. For ee and e/x 
spectra, the signal is bigger than the background. For yU/i spectra, the signal 
is smaller than the background, but with the help of like-sign subtraction. 
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it is possible to also observe the /x/i signal from open charm decay. In chap- 
ter we show the signal and background for the single inclusive electron 
and muon spectrum. The signal from single inclusive leptons reflects the 
shadowing effects less effectively, and the signal-to-background ratio there is 
much smaller than that for the lepton pair case. Finally in chapter |4.12| we 
give the discussion and summary. 



4.2 A Recent Attempt to Extract Nuclear Gluon Shad- 
owing 



Recently there appeared an attempt ||174|| to extract nuclear gluon shadowing 
from the high statistics data on Sn and C taken by NMC | 175 |. 



Parton structure functions depend on Bjorken x, the momentum fraction 
of the parton in the frame where the nucleon (or nucleus) is moving very 
fast, and Q^, the momentum scale at which the parton structure function is 
being resolved. The evolution equations for the gluon number density, 
G{x), and the quark number densities, qi{x) {i from 1 to 2n/), are given by 
the Altarelli-Parisi equations ||176|| as 



dG{x) 


a 


/•I dy 


dt 


27r 


Jx y 


dqi{x) 


a 


r^dy 


dt 


277 


L y 



Y,PG,,{z)qM + PGG{z)G{y) 



(4.2) 



where z = x/y, and t = \n{Q^/Ql). There are four splitting functions, 
PGg^iz),PGGiz),Pgg{z),Pg^g^{z) 3.11(1 Pg^G{z), Hud a function Pab{z) is basi- 
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cally proportional to the probability that a parton B splits to a parton A 
with fraction z of the momentum of the parent parton B. 

Based on the Altarelli-Parisi equation ( ^4.2| ) for quarks, the evolution 
equation for the quark distribution 



(4.3) 



is therefore given as 



dt 



a dy 
27r Jx y 



a dz 
2tt Jx z 



J2e^P,,,^{z)q,{y) +J2e^P,^G{z)G{y) 



^ X ^ 



a 
2^ 



(^Pg,*F2 + xJ2elPgG*G'^ 



(4.4) 
(4.5) 
(4.6) 



The above two convolutions are defined as 



P,G*G 



Pqq * F2 



f -P,g{z)G{-) = r ^[z' + (1 - zm^) (4.7) 

Jx Z Z Jx ^Z Z 

C dzP,,{z)F2{-) , 

J X Z 



•1 dz. 



1 N'^ - I 
dz- " 



2Nr 



1 + '^' 



,x. 



4 fi 



3 Jx 



dz 



il + z')F2{x/z)-2F2{x) 3 _ X. 

(l-z) 2 ^ ^ ^^2' 



,(4.8) 



with the divergence at ^ = 1 being regularized by the definition 



(1-^)+ Jo (l-z) 



(4.9) 
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Using the approximation 

f'dz[z^ + il-zf]-G{-) ~ ^G{2x) , (4.10) 

J X Z Z O 

one can relate the second term in eq. (|4.6|) , Pgc * G, directly to the gluon 
density function G{x). Note that the first term in eg. ( |4.6|) , Pgg * F2, is only 
related to F2{x), which can be computed from data. 

Thus from the dependence of the quark distribution F2, one can ex- 
tract some information on the gluon density G. In the above, one neglects 
the effects of higher order corrections and the non-linear recombination terms 



131 1 . The non-linear term probably does not have a large magnitude unless 



the parton momentum fraction x is extremely small |[L31| , |166| , |177 



In Figure [4.1| ||174|| , the so-extracted ratio of gluon density of tin to carbon 



r(x, Q^) = G^"{x,Q^)/G'-^{x,Q'^) is plotted against the gluon momentum 
fraction x in small x range 0.01 < x < 0.1, together with the F2{x, Q"^) data 
from which the ratio r{x, Q"^) is extracted. The ratio of gluon density shows 
both the shadowing effect in the lower range of x and the ant i- shadowing 
effect in the higher range of x, similar to the case of quark shadowing. 



4.3 The Evolution of Nuclear Parton Distributions 

There are mainly two steps to determine nuclear shadowing effects. First 
we need to know the parton distribution functions in nuclei at some initial 
momentum scale Qq. Given the input nuclear parton distributions at the 
scale Qo, "we then calculate the parton structure functions at any values 
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Figure 4.1: The ratio r(x) = G^"' (x) / (x) of tin to carbon gluon density, 
together with the ratio of structure function, = F2^{x) / F2 { 



function of x, from [174 . 
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of scale Q according to evolution equations. As long as Q is within the 
perturbative QCD domain, evolution equations are known from pQCD 



176, 131 



The input parton structure functions, G{x,Ql) and qi{x,Ql), contain all 
non-perturbative nuclear effects at the scale Qo- They are thus, in principle, 
not calculable within pQCD. However, it is not hard to imagine that nuclear 



shadowing effects occur to nuclear parton distributions ||178|| . We can imagine 
that shadowing effects occur when partons from different nucleons in nuclei 
begin to overlap spatially. For a nucleus with momentum p per nucleon, 
the longitudinal size of the parton with momentum xp is ~ l/{xp). The 
longitudinal separation of different nucleons is ~ 2rom/p, where m is the 
nucleon rest mass, and tq is the nucleon radius. Thus shadowing effects 
occur when x < l/(2rom). Shadowing effects grow as partons overlap more, 
until the partons are overlapping completely. This saturation of shadowing 
effects happens when the longitudinal size of the parton is about equal to the 
longitudinal size of the nucleus, i.e. x ~ l/(2i?^m), where Ra is the radius 
of the nucleus. 

evolution equations of parton density functions are usually given by 
Altarelli-Parisi equations (AP) ||176|| , as expressed in eqs. ([4.2|). They can 
be rewritten separately for the gluon momentum distribution xG, the flavor 
singlet sea quark part xS, and the flavor nonsinglet valence quark part xV: 

+zPGG{z)yG{y,Q'')\ , (4.11) 
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zPgg{z)yS{y,Q'^ 



am f'dy 
2iT Jx y 

+zPgG{z)yGiy,Q'] 



d 



■xV{x,Q^ 



['^zPggiz)yViy,Q' 
Jx y 



dQ^^' ' 27r Jx y 

where the sea quark momentum distribution xS = xUs,xds or xsg 
The four sphtting functions are given as ||176 



P,g{z) 
Pgg{z) 



1 - z 



'1-z) 



z(l-z) 



p,,i^) 



Pg,{z) 



A 1+z^ 



3 (1 - , 
41 + (1-^)2 



+ 2(5(1 - z) 



(4.12) 
(4.13) 



(4.14) 



(4.15) 
(4.16) 
(4.17) 



3 z 

where {1 — z)^^ is defined in eq. (^^). 

The AP-evolution means that as increases, we can see more small-x 
partons and less high-x ones, because high-momentum partons lose energy 
by radiation. Thus the parton density can become very large at sufficiently 
small X. However, as the parton density becomes huge at small x, the AP- 



evolution equations have to be modified | ]131| ] because gluons then have more 
probability to recombine to reduce the numbers of partons in the nuclei. 

modified the AP-equations by calculating the 



Mueller and Qiu [TTl 



recombination probabilities for gluons to go into gluons or into quarks. The 
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AP-evolution equation( [4.11| ) for gluons and equation ( [4 .121 ) for sea quarks are 
modified to the following AP+MQ-equations: 



+zPGG{z)yG{y,Q': 



(4.18) 



d 



zP,g{z)yS{y,Q') 



am rdy 
2tt Jx y 

+zPgG{z)yG{y,Q^) 

n 

- 

where the nuclear number density n ~ l/(7ri?^), and the contributions to the 



i^^L^M') + RfUM')] ,(4.19) 



evolution from recombinations of gluons into gluons or into quarks are 



- 1 



IT 



'dy r 

X y 



yG{y,Q') 



xG{x, 



2a{Q^)Tf f^dy 



TT 



X y 



-Z-fFG{z)yGHT{y,Q'^ 



(4.20) 
(4.21) 
(4.22) 



In the above Tj = 1/2, Ca = N = 3, ^fg{z) = -2x + 15x^ - 30x^ + 10x\ 



(4.23) 



gives the evolution for the higher- dimensional gluon distribution Ght- This 
distribution function is basically not calculable [|178|| . Since it comes from 
two-gluon interaction, it has the form, at some scale Qq, as 



yGHT{y,Ql) = KHT yG{y,Ql] 



(4.24) 
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where Kht is a constant. 

From the above AP+MQ-equations, the modification terms from two- 
gluon recombinations to the AP-equations are basically of the form 

a\yGf/{R^Q') . (4.25) 

We take gluons for example, and roughly represent the modified evolution 
equation (|4 .181) as 

Q ~ a—xG{x,Q^)-a'-^[xG{x,Q^ . (4.26) 

The modification term, the second on the right-hand side in eq. (|4.26|) , there- 
fore slows down the Q"^ evolution of the gluon density function. To make a 
rough estimate of the onset of non-linear recombination effects, we equal the 
two terms in eq. (|4.26|) , the usual AP-evolution term and the gluon recombi- 
nation term, and get the condition as 

^G{x, Q^) ~ ~ ^ . (4.27) 



47r2 



na Tca 



The modification will be significant only when x is so small that the above 
condition is reached. If we take R = Ifm, Q 2GeV, the above condition 
will be reached only for extremely small x where xG{x, Q"^) ~ 100. Thus the 
recombination is likely to be a very small correction to the usual evolution 
equations of parton structure functions. 

Modified evolution equations (AP+MQ) was first numerically solved by 
Qiu ||178|| , and then by Eskola ||166|| . Eskola assumes a shadowing ansatz at 
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an initial scale Qo = 2GeV for heavy nuclei with A ~ 200, then evolves shad- 
owed Duke-Owens parton distributions ||144|| according to both AP equations 



and modified evolution equations (AP+MQ). The MQ-modification slows 
down the increase of the parton densities with Q^, and therefore shadow- 
ing effects go away slower with the increase of than in the AP-evolution 
case. However, the difference between the AP-evolution and the AP+MQ- 
evolution with this particular choice of input parton distributions is found 
to be only on the order of 10% |166|, thus the mainline of the evolution of 



parton structure functions is given by the usual AP-equations. 

4.4 Nuclear Shadowing in Open Charm Production 

We use two common shadowing scenarios: the parameterization which does 
not depend on scale or component (HIJING's scaling) []70|, and the scale- 



dependent, component- dependent shadowing (Eskola's shadowing) ||166|| . The 
production for cc in one A + B collision at a fixed impact-parameter b is given 
by H 



J 2 J J = ^ / d'^ry2^a'ra/Aixa,Q^,r}XbTh/Bixb,Q'^,b-r}—^ ,(4.28) 
dp\dyidy2 ■'7^ dt 

where 

^aM(x,Q^F) = TA{f)fa/N{x,Q'')Ra/A{x,Q'',r} (4.29) 

is the nuclear parton density function in terms of the known nucleon par- 
ton structure functions fa/N{x,Q'^), the nuclear thickness function TA^r) = 
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/ dznA{\J + f^), and the unknown impact-parameter-dependent shadowing 
function Q^, t). The conventional kinematic variables Xa^h are the in- 

coming parton light cone momentum fractions, t = —{pb—pi)'^, and the final 
c, c have rapidities yi,y2 and transverse momenta pi± = — P2± = Ri- Next- 
to- leading order corrections to the lowest order parton cross sections dcTab/dt 
are approximately taken into account by a constant K factor ||151| , |152|| . We 



use the recent MRSA [ p.45|| parton structure functions. From a fit to low 
energy open charm production data in pp collisions, we fix rric = 1.4 GeV, 
K = 3, = s/2 for gg cc and s for qq cc. This choice of parameters 
leads to a charm pair cross section of 340/i6 for pp collisions at RHIC. For 
impact-parameter-averaged collisions, the integral over the transverse vector 



in eg. ( |4.28| ) leads to a factor: 

BAR,/Bixb,Q^)Ra/Aixa,Q^)/af/ , (4.30) 

where is the inelastic B + A cross section and Ra/A{xa, Q^) is the impact- 
parameter-averaged shadowing function. 

Therefore for cc production in impact-parameter-averaged p + Au colli- 
sions, we have 

dp^ dyidy2 " ^^^^''f''/^^^''^Q'^'>^^f^/^^^''^Q'^'>^btAui^b,Q'^)-^ 



where Rb/Aa{xh, Q^) is the impact-parameter-averaged shadowing function 
for the parton b in Au nucleus. 

The two different shadowing functions are plotted in Figure [4.5|. It is clear 



that HIJING's scenario has much stronger shadowing than Eskola's scenario 
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in almost all regions of interest, and Eskola's scenario has anti-sliadowing in 
X > 0.1 region. These features are reflected in the final e/x spectrum, as we 
will show in the following chapter. 



4.5 e/i Spectrum from D/D Decay 

Since the mass difference of e and yU has little effect on the semileptonic decay 
of charmed meson, the spectra of dilepton pair ee and /i/i from charm decay 
have shapes almost identical to the e/x spectrum. Since we include both e'^n' 
and e^/i+ pairs, and assume the average branching ratio for — > /i + X is 
the same as for D e + X, the e/i signal is twice as large as the ee and /i/x 
signal. 

To compute the DD pair distribution, a hadronization scheme for c ^ D 
must be adopted. In e+e^ — cc ^ DDX hadronization can be modelled 
via string fragmentation or fitted, for example, via the Peterson fragmenta- 



tion function (see ||141|| ). The final D carries typically only a fraction ~ 0.7 
of the original c momentum. However, in pp collisions, charm hadroniza- 
tion is complicated by the high density of partons produced during beam 
jet fragmentation. In this system recombination or coalescence of the heavy 
quark with co-moving partons provides another mechanism which in fact 
seems to be dominant at least at present energies. The inclusive pp DX 
data is best reproduced |p.41| , |179|| with the delta function fragmentation 



D{z) = 5(1 — z) in all observed Xf and moderate pt regions. This obser- 
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vation can be understood in terms of a coalescence model if the coalescence 
radius is Pc ~ 400 MeV. We assume that hard fragmentation continues to 
be the dominant mechanism at RHIC energies, and therefore no additional 
A-dependent effects from hadronization will arise. However, this assumption 
must be tested experimentally. Therefore our main dynamical assumption 
is that the pp D transverse momentum distributions can be accurately 
reproduced from the QCD level rates assuming hard fragmentation, as at 
present energies ||179|| . At RHIC this can be checked either via single inclu- 



sive leptons 1 39 1 or directly via Ktt | 180 |. 

With this assumption, the impact-parameter-averaged DD pair distribu- 
tion in p + A is given by 

dNP^ KA EsE^ 



dpldysdy^ (jPf E1E2 

>^J2^bfb/NiXb,Q'^)Xafa/NiXa,Q'^)Ra/AiXa,Q'^)^^^ .(4.32) 
b,a ^'^ 

Note that there is no kick here, and D and Z) in a pair have opposite 
azimuthal directions. 

After the above spectrum dN/dp^dy^dy^ for the charmed mesons is cal- 
culated, we use Monte Carlo programs to simulate their semileptonic decays 



181|| . We generate D/D events according to the above spectrum. For each 



event, we decay one charmed meson in its rest frame to an electron, and 
decay the other charmed meson in its rest frame to a muon isotropically, ac- 
cording to the parameterized semileptonic energy spectrum. We then boost 
the electron and muon back according to the momentum of its meson parent. 
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We then calculate the invariant mass 



M = ^{Pe+V,? (4.33) 

and the rapidity 

y = tanh-i[(pll + E^)] (4.34) 

for the pair and have the final e/i pair spectrum dN/ dMdy. 

In Figure ^3(a) we plot the absolute rate da^^ /dMdy from the above 
Monte Carlo calculation, and it includes both e^fi^ and e~fi~^ signals. The 
solid, dashed and dot-dashed curves represent the non-shadowing, Eskola's 
shadowing, and HIJING's shadowing case, respectively. From top to bottom, 
the three sets of solid, dashed and dot-dashed curves correspond to different 
pair-masses M = 1, 2, 4 GeV. In Figure ^^b), the ratio of the curves in (a) 
from two shadowing cases to that from the non-shadowing case is plotted 
at pair- masses M = 1 (solid), 2 (dashed), 4 (dot-dashed) GeV. The up- 
per (lower) set of three curves corresponds to Eskola(IIIJING)'s shadowing. 
The curves in Figure |4.2| (b) resemble the shapes of shadowing curves in Fig- 



ure [4.5| . Note that in Figure |4.5| the abscissa is in logarithmic scale and in 
Figure |4.2| the abscissa is in linear scale. We see that the shadowing effect 
is much stronger in positive pair-rapidity region. This is expected since in 
p + A collisions positive pair-rapidity corresponds to smaller x from nucleus 
A and thus deeper shadowing. For Eskola's scenario some ant i- shadowing 
effect shows up in the negative pair-rapidity region. One may notice that in 
Figure [4.2K b) we have a narrower rapidity region at higher invariant masses. 
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That is because at the edges of the rapidity region we have fewer Monte Carlo 
events, and we discard the edges to restrain big statistical fluctuations. 

4.6 Scaling 

Since gluon fusion dominates, there is an approximate scaling in terms of the 
light cone fraction xa that one of the gluon carries from nucleus A: 

InxA = -Vcc + ln(Meg/yi) • (4.35) 

While the lepton pair-mass M and rapidity y fluctuate around a mean 
value for any fixed M^g and ycc, on the average 

(y) ~ (yc) , (4.36) 

and (Mcc) can be well approximated by a linear relation 

{M,,)^(3M + Mo (4.37) 

over the lepton pair-mass range 1 < M < 4GeV. For the delta function 
fragmentation case, (3 ~ 1.5 and Mq ~ 3.0 GeV are determined by the D- 
meson decay kinematics. In Figure [4.3| we plot the average M^g as a function 
of M, the mass of the correlated lepton pair from the cc decay with the delta 
function fragmentation. The above relation between the average masses is 
determined from the following histogram plot generated from Monte Carlo 
calculation, which showed the event number density on the plane of charm 
pair-mass Mcc and the mass M of the dielectron from the charm semileptonic 
decays (the unit of mass is GeV in the figure) . 
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Figure 4.2: (a) The production of opposite-sign e/i pairs as a function of 
pair-rapidity y at pair- masses M = 1,2,4 GeV; (b) The ratio of shadowing 
curves over non-shadowing curves in (a). 




Figure 4.3: The average mass Mcc as a function of dilepton mass M. 
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dN/d(Mcc)d(Mee) 



Figure 4.4: Probability distribution of M, the dilepton mass from charm 
pair decay, as a function of Mcc in p + Au coUision. 
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Therefore, 

InxA ~ -l/ + ln[(/3M + Mo)/v^ . (4.38) 

The ratio of the dilepton dN/ dMdy spectra in pA to scaled pp for different 
pair-masses is thus expected to scale approximately as 

<^ (M, y = -\nxA + ln[(/3M + Mo)/v^ ) 



1 dN'pJ 



Rg/AixA,Q'^iPM+MoY/2) , (4.39) 



where u ^ AaZlol^ ~ A'l\ 

In order to test gluon dominance and the accuracy of the above approx- 
imate scaling, we compare in Figure |4.5| the gluon shadowing function to 
the above dilepton ratio from the Monte Carlo calculation. The upper solid 
curve is Eskola's gluon shadowing | |166| | for = 10 GeV^, and the lower 



solid curve is shadowing from HIJING |7^. The other six curves are ratios 
of dilepton dN/ dMdy spectra of shadowed p + Au over those from unshad- 
owed p + Au, as given by eq.( [4.39| ) as a function of the scaling variable xa 



for three different dilepton masses. In Figure ^3] we first plot all the ratio 
curves in terms of reversed pair-rapidity —y, then we shift the ratio curves 
at M = 1, 2, 4 GeV to the left by 3.79, 3.51, and 3.10 respectively. 

Overall, with an approximate scaling the ratio of spectra from p + A to 
those from pp reflects the nuclear shadowing function well. We conclude 
that this ratio can therefore serve to map out gluon shadowing in nuclei. 
Note that in the case of Eskola's shadowing, even the dependence of the 
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Figure 4.5: Dilepton dN/dMdy ratio curves (shadowed over unshadowed) 
as a function of the scahng variable logioa;^ at different masses, and the 
comparison with shadowing curves. 
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shadowing function is visible through the rise of the ratio curves with M in 
the small x region. 
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4.7 Original Backgrounds 



It is important to estimate the dilepton background to see if the proposed 
signal is experimentally feasible. In this chapter we estimate some dominant 
backgrounds without taking any suppression or kinematical cuts. In chap- 
ter [4.8| we calculate the signal and background dileptons for the PHENIX 
detector P^, taking into account the detector geometry and specific kine- 



matical cuts. 

The electron background comes from Dalitz decay of 7r° mainly, and also 
from Dalitz decay of higher mass particles and resonances. Also we have 
electron background from photon conversion. For dielectron measurements, 
7r° Dalitz decay background (vr^ e^e~7) is quite problematic. The 

correlated pair from Dalitz decay has an invariant mass smaller than the pion 
mass, so ideally we can eliminate all these electrons and positrons by invariant 
mass constraint. However, electrons with small momenta may be easily lost 
in the strong magnetic field in the detector, in addition, we will not have An 
solid angle coverage. Although we can reject many of these leptons by using 
techniques like small-angle cuts for dielectron pairs, we still have background 
from 7r° Dalitz decay, and this is also a problem for e/x measurements. The 
main muon background comes from the random decays of charged pions and 
kaons. They can decay in the free space before the hadron absorber, and they 
can also decay in the absorber or leak through. According to the momentum 
resolution in the PHENIX muon detector [aiOp) = O.4°,a(0p) = 1.6°) ||, 
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it is very difficult to differentiate these muons from muons coming directly 
from charm decay. 

Let us estimate the dielectron numbers in a central p + Au collision at 
200 GeV/A at RHIC. Given N as the number of vr^s, P as the branching 
ratio of 7r° —>■ e+e~7, the average number of combinatoric (uncorrelated) 
dielectron pairs is then A^(A^ — 1)P^. The HIJING Monte Carlo calculation 
| 70| gives about 38 vr^s, and the branching ratio of 7r° Dalitz decay is about 
1.2%. Therefore the average number of dielectrons originating from the 7r° 
background is 0.2. Compared to the next largest source (dileptons from the 
D/D decay), this background is enormous. Using HIJING's shadowing, we 
have 0.05 D/D pairs. With the average branching ratio of D — > e + X at 
about 12%, the dielectron yield from charm decay is about 0.0007, a factor 
of 300 below the original background from vr*^ Dalitz decay. For dimuon 
numbers in a central p + Au collision at RHIC, HIJING gives about 34 vr+s 
and 7r~s each, and they almost all decay to muons. Thus, the original dimuon 
background from charged pions is 1.1 x 10^, which is enormous compared with 
the dimuon yield from charm decay (also about 0.0007). For opposite-sign 
e/i numbers, the original background originating from vr" and 7r+(7r^) is 31., 
compared to e/i from charm decay (about 0.0014). 

The above provides us with a rough idea of the opposite-sign lepton pair 
backgrounds before we suppress them. We must suppress them in order to 
observe signals from charm decay and thus probe nuclear shadowing. 
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4.8 Backgrounds Entering the Detector 



In this chapter we study the backgrounds entering the detector by includ- 
ing detector-related effects: the suppression factor for leptons, proposed 
PHENIX detector geometry, and a specific kinematical cut for the events. 
These are based on the RHIC report . Then we compare the backgrounds 



to the lepton pairs signal from charm decay. 

The electrons from vr'^ Dalitz decay could be suppressed by the small- 
angle cut for dielectron pairs since they have an invariant mass less than 
the pion mass. We suppose we could reject 90% of the 7r° Dalitz electrons 
|18!2|. We also have electron background from photon conversions, which has 



a softer spectrum. We suppose we have the same rejection rate of 90% for 
the conversion electrons. Each of the two photons decayed from vr" has a 
0.5% probability of conversion to electron/positron pair, and the 7r° Dalitz 
decay branching ratio is 1.2%; as an upper limit we will simply double the 
electron background from 7r° Dalitz decay to include the contribution from 
photon conversion. From HIJING program, in terms of number of electrons 
from tt" Dalitz decay, there is approximately 40% more electrons from K~ 
decay. We also assume a rejection rate of 90% for these electrons based on 
the estimates of free space decay probability and angular resolution for the 
vertex. Thus, the number of background electrons after rejection is 24% of 
the number of ti^ Dalitz electrons, which means we have a suppression factor 
of about 4 for the original background electron number from 7r° Dalitz decay. 
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We generate the background electrons from the electron spectrum calculated 
from HIJING program (using proper weights for Dahtz and K~ Dalitz 
events), since electrons from K~ have a harder spectrum, and thus a greater 
probability of passing the kinematical cut. 

The background muons can come from decays in the free space before 
the hadron absorber, which is about 32 centimeters(L) from the nominal 
vertex. For a charged pion with mass m, energy E, and proper decay time 
T, the probabihty of the decay in free space is about Lm/rE. We would 
use an energy cut of 2 GeV for muons, therefore this probability is less than 
0.3%. There is also a background from particles that have leaked through 
the absorber, which decay to muons in flight inside the muon arm. Prom the 
RHIC report, for charged hadrons above 1 GeV in the muon arm for a central 
AuAu collision at 200 GeV, the leakage number is 4.8, which is about 0.1% of 
the total number of charged pions. Their decay to muons is also suppressed 
by the long proper lifetime r and Lorentz 7 factor as above, and those which 
have not decayed will mostly be differentiated by the muon identifier. So 
the probability of finding a muon above 1 GeV from the decays of leakage 
particles is about 0.003%. Considering the factor from the geometry and 
energy cuts explained below, the probability from leakage before these cuts 
is also 0.3%. The charged hadrons could also decay in the absorber during 
the shower. After all, we take 1% of the muon number from charged pion 
decay as the estimate of the muon background passing the hadron absorber 
and muon identifier, which means we have a suppression factor of 100 for the 
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original background muon number. We generate background muons from 
the muon spectrum calculated from HIJING program, which is mainly from 
TT"*", 7r~ two-body decay, and, in very small amounts, from decay of heavier 
particles such as K'^ and K". 

For the PHENIX detector geometry, the electron arm barrel covers elec- 
tron pseudo-rapidity range 

- 0.35 <rie< 0.35 , (4.40) 

and azimuthal angle range 

0e e ±(22.5°, 112.5°) . (4.41) 

The muon arm endcap covers the polar angle range 10° < ^ < 35° and almost 
the full azimuthal angle. Thus, we only require muon pseudo-rapidity range 
to be 

1.15 < 77^ < 2.44 . (4.42) 

For the kinematical cut used to improve the signal-to-background ratio, we 
take 

Ee > IGeV , (4.43) 
> 2GeV . (4.44) 

We also require the relative azimuthal angle of the lepton pair to be 4>i+i- > 
90°, and this improves the ratio by a factor of 2. We find the following 

percentages for the numbers of leptons which pass the above energy and 
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pseudo-rapidity requirements: 1.1% for electrons from charm decay, 2.7% 
for muons from charm decay, 0.02% for background electrons, and 1.0% for 
background muons. The small percentage for background electrons is due to 
the small mass of the particles which make Dalitz decays, and this greatly 
lowers the electron background. Recently, a second muon arm at the opposite 
end is proposed to increase the muon rapidity coverage. 

Considering the above suppression factors resulting from the geometry 
and the kinematical cuts, we calculate the signal and the backgrounds enter- 
ing the detector. First, the ee, e/i and /i/x signals entering the detector from 
charm decay are calculated in three different shadowing cases. The dN/dM 



plots of different shadowings are shown in Figure [4.6| , and the dN/dy plots 



are shown in Figure |4.7| . In Figure |4.6| , results from non-shadowing case 
(solid), Eskola's shadowing (dashed), and HIJING's shadowing (dot-dashed) 
are plotted. Figure (a), (b), (c) are for opposite-sign ee, e/x and /i/i, respec- 
tively. In Figure [4.7| , the left, middle and right set of three curves correspond 
to opposite-sign ee, e/x and /i/i signals. Note that the e/i and /i/x yields are 
scaled down by a factor of 4 and 15, respectively. Due to detector geome- 
try and kinematical cuts, ee, e/x and /i/i spectra cover pair-rapidity regions 
centered at about 0, 1 and 2, respectively (the second PHENIX muon arm 
can reach pair-rapidity region centered around —1 and —2). As shown in 
Figure ^]2|(b) for the signal without cuts, the shadowing effect is stronger in 
larger pair-rapidity region. 

We then generate background electrons and muons based on HIJING cal- 
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Figure 4.6: dN/dM spectra for lepton pairs entering the PHENIX detector 
for different shadowing cases. 




Figure 4.7: dN/dy spectra for lepton pairs entering the detector for different 
shadowing cases. 
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culation, and calculate the ee, e/i and jjjj, backgrounds entering the detector. 
We compare the backgrounds with the calculated signal from the HIJING 
shadowing case. The dN/dM plots of signal and backgrounds are shown in 



Figure |4.8| , and the dN/dy plots are shown in Figure In Figure |4.8| , 
the solid curves represent the signals from open Charm decay — labeled CC 
signal. The label C refers to the lepton originating from open Charm decay, 
label D refers to the electron originating from Dalitz and photon conversion, 
and label R refers to the muon originating from Random decay of pions and 
kaons. In (a), the dielectron backgrounds are represented by the dashed curve 
DC and the dot-dashed curve DD; and both are scaled up by a factor of 100. 
In (b), the opposite-sign e/i backgrounds are represented by the dotted curve 
CR, the dashed curve DR, and the dot-dashed curve DC. In (c), the dimuon 
backgrounds are represented by the dashed curve RR and the dot-dashed 



curve RC. The notations in Figure [4.9| are the same as in Figure ^^8 . 

We find that the signal-to-background ratio for ee is very large, thus the 
dielectron signal from open charm decay is the easiest to extract. That ratio 
falls to about 2.5 for e/x, and about 1/4 for yU/i. Note that we have not 
included the backgrounds from Drell-Yan and thermal production that will 
most change the above signal-to-background ratio for ee. 

One could further suppress the lepton pair backgrounds by like-sign sub- 
traction, especially in the /z/x channel. Since the particle production has a 
fairly good charge symmetry, the background opposite-sign lepton pairs are 
almost the same as the like-sign pairs. The signals from open charm de- 
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Figure 4.8: Charm signal and backgrounds entering the detector as a func- 
tion of the lepton pair-mass M. 
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Figure 4.9: Charm signal and backgrounds entering the detector as a func- 
tion of the lepton pair-rapidity y. 
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cay are almost all opposite-sign pairs, so we could subtract like-sign pairs 
from opposite-sign pairs to suppress the backgrounds. The efficiency of the 
like-sign subtraction depends on the charge and isospin asymmetry of the 
collision. With the help of like-sign subtraction, the e/i and fifi signal could 
be observed. 



4.9 Single Inclusive Lepton Spectrum 

The kinematics here are considerably less restricted, so the relation between 
the kinematic variable E, rj for the single lepton and the gluon Bjorken vari- 
able xa are more indirect. Thus, the single inclusive lepton spectrum reflects 
the nuclear shadowing effects less effectively. For the signal from open charm 
decay, the ratios of electron pseudo-rapidity {rj = l/21n[(l -|- cos6')/(l — 
cos^^)]) distributions dN/drj from the two shadowing cases over those from 
the non-shadowing case are shown in Figure [4.10| . The solid(dashed) curve 
represents HIJING(Eskola)'s shadowing. 

Also the signal-to-background ratio in the detector is not as good as that 
for the lepton pair case. The estimated ratio is approximately 3 for single 
electrons and 1/22 for single muons (with HIJING's shadowing), compared 
with about 260 for ee, 2.5 for e/x and 1/4 for /i/i. Given the above small 
signal-to-background ratio, it is difficult to observe single muons. Since the 
single electron and muon spectrum cover pair-rapidity ranges around and 
0, 2 respectively, we can only observe mid-rapidity single electrons. While 
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Figure 4.10: Ratios of single electron signals from the two shadowing cases 
over those from the non-shadowing case as a function of the pseudo-rapidity 
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for the the three pair spectra, they cover rapidities around 0, 1 and 2, and 
have reasonable signal-to-background ratio. So we conclude that although 
the single electron spectrum signal could be seen as an easy initial measure 
of shadowing effects, the shadowing effects will be seen much more clearly in 
the lepton pair spectrum than in the single spectrum. 



4.10 Cronin Effect and Energy Loss 

Cronin effect and energy loss will weaken the back-to-back type of correlation 
between the two charm mesons in the pair, and thus weaken the scaling 
behavior among the dilepton dN/dMdy ratio curves showed in Figure ^75. 



However, in the following they are estimated to lead to distortions of up to 
10% m. p + A collisions at RHIC energies. 

From studies of the nuclear dependence of transverse momenta of J/%1) 
production, the typical increment, 5p\{A), due to multiple collisions is lim- 
ited to ~ 0.34 GeV^ even for the heaviest nuclei ||183|| . This momentum 



spread is shared between the c and c, and distributed approximately as 

g{6p^) = e-^^'i/^VTTA' , (4.45) 

where ~ O.lTGeV^. For energy loss, we assume that the charm quark 
loses an energy 5E in the lab frame, and the charm quark is reduced by 
(1 — e), where 



e = (5E/[m_L sinh(?/ + yo)] 



(4.46) 
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Combining these two effects, we may write for tlie cliarm quark tfiat the final 
transverse momentum is 

P± = (pi"' + ^P±)(l-e) . (4.47) 

Assuming that the D-meson spectrum without these two effects can be ex- 
pressed as 

F(p^) = (fN/dm\ , (4.48) 

then it becomes 

F'{p±) = I F{pT)g{5p^)d5pi. . (4.49) 

with 

F(pi_) (X e-"^^ ; a ~ 1.3 GeV^ (4.50) 

from HIJING. Expanding the convolution to lowest order, the relative change 
of the D-meson spectrum is given by 

F'(p±) /F{p^) ~ 1 - ap^e + a" A74 , (4.51) 

This is also the relative change of the cc pair spectrum when M^c — 2m± in 
case of yi = y2, P±i = P±2- Take 

mx~3GeV , (4.52) 

5E ^lOGeV , (4.53) 

coshyo ^ 100 , (4.54) 
then the relative change in the pair spectrum is estimated to be 

F'{p±)/F{p±) ~ 1 - 0.1(e-^ - 1) . (4.55) 
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4.11 Effect of Energy Loss in ^4 + ^4 Collisions 

In heavy ion experiments, dileptons from a few GeV region were suggested 
to be possible signals of the quark-gluon plasma; however, they were later 
claimed to be dominated by the dileptons from open charm decay. Nev- 



ertheless, it has been suggested recently |^9[ that energy loss is an important 
factor which could soften greatly the dilepton spectrum from open charm 
decays. 

In y4 + y4 collisions, open charm must pass through the dense parton 
medium formed after the initial collision; it will therefore lose a certain 
amount of energy, typically 2GeV/fm. Moreover, the two charm mesons 
in a pair tend to have opposite azimuthal angles, so the probability that 
both of them can escape without substantial energy loss is very small. Con- 
sidering that we are interested in dileptons with energies of a few GeV, the 
dilepton spectrum from correlated open charm (also open bottom) decay 
will be greatly softened compared to the conclusion in and is possibly 
comparable to the signal from the quark-gluon plasma. 

In order to get a better understanding of the magnitude of the softening 
of the charm dilepton signal, we can take into account the following facts: 
in a finite size nucleus, partons lose energy in finite number of collisions, not 
by the simple linear energy loss formula. Therefore the energy loss has more 
fluctuations than in the linear assumption case. Secondly, the dilepton p± 
spectrum depends crucially on the angular distribution, thus it is important 



151 



to calculate the open charm production to the next-to-leading order in order 
to take into account the fact that the open charm pair are not always back- 
to-back azimuthally. 

In the mid-rapidity region in p+A collisions, there are much fewer partons 
than in the A -|- A case, and the most energy loss occurs in the target rapidity 
region where y << 0. This energy loss effect was estimated in the above 



chapter 4.10 and shown to be very small. 



4.12 Discussion and Summary 

We use Delta function as the charm quark fragmentation function to D me- 
son. The above analysis depends on the validity of this basic assumption 
regarding the hard fragmentation of charm quarks. Although this assump- 
tion is being supported by the low energy pp data ||141| , |179|| , it should be 



checked explicitly via the single inclusive measurements of D production at 
RHIC. 

In this study, we included only the geometrical part of the experimen- 
tal acceptance in our calculations. In order to get the detailed background 
spectra for photon conversions, leakage muons, etc, the full detector simu- 
lation should be carried out. Since the geometry is asymmetric for p + Au 
collisions, energy loss will shift the rapidity of the particles toward the target 
rapidity region, which is an additional shift to the larger rapidity shift due to 
shadowing effects. Jet quenching could also change the rapidity distribution 
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slightly if the nuclear density is high, particularly for the high p± part. We 
must separate contributions due to other effects such as energy loss in order 
to clearly observe the interesting shadowing curves shown in Figure 



We choose p + Au since it has much smaller combinatoric background. 
It also demonstrates the shadowing effects more clearly in the dN/dy{M) 
spectrum for the lepton pairs as shown in Figure 



In summary, we calculated the lepton pair spectra from open charm decay 
in two different shadowing scenarios. By scaling the ratios for different mass 
ranges according to eq.(^]3^), we showed that the dilepton rapidity depen- 
dence of those ratios on xa reproduces well the underlying gluon shadowing 
function defined in eq. ([4.1|) . Finally we showed that the measurements re- 
quired to extract the gluon shadowing are experimentally feasible at RHIC. 
The signal can be seen in p + Au collisions, some even before using like-sign 
subtraction. 

We conclude by emphasizing the importance of determining gluon shad- 
owing in p + y4 to fix theoretically the initial conditions in A + A. In A + A 
the open charm decay is regarded as an annoyingly large background that 
must be subtracted to uncover the thermal signal. In p + A that background 
becomes the signal needed to determine the incident gluon flux in A + A. 
The continuum charm dileptons in p + A at RHIC are likely to provide a 
unique source of information on the low xa nuclear gluon structure at least 
until HERA is capable of accelerating heavy nuclei. 



153 



Chapter 5 
Outlook 



In this thesis work, we have investigated the open charm production in nu- 
clear coUisions and whether it can serve as a probe of the quark-gluon plasma 
phase. We found that the additional open charm production from the pre- 
equilibrium stage is negligible if we allow strong correlations between the 
space-time coordinate and the momentum of minijet partons. Contrary to 
an earlier result, which claimed an enormous enhancement in the open charm 
production given the formation of the quark-gluon plasma, we showed that 
open charm is not likely to be a good probe of the plasma phase formed in 
relativistic heavy ion collisions. 

We then proposed to use the open charm as a probe of the initial parton 
structure function in nuclei, and thus turned our attention from probing the 
quark-gluon plasma to probing the nuclear shadowing effect. The dominance 
of the open charm in lepton-pair spectra with a few GeV energy makes it 
possible to observe the open charm via high- mass lepton-pairs in p + A colli- 
sions. To make a better connection with ongoing experimental projects, we 
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then estimated the signal-to-background ratio for the proposed lepton-pair 
measurements for the PHENIX detector at RHIC, and showed that these 
measurements are feasible. 

In our study, we find that the pre-equilibrium charm production depends 
strongly on the space-time and momentum correlation of minijet partons. 
Our model is a minimal modification of the Bjorken correlation by taking into 
account the uncertainty principle. The exact strength of the correlation of 
minijet partons is yet an open question, and needs to be further investigated. 

The assumption on the hard fragmentation function of charm quarks 
needs to be tested in high energy hadroproductions. The latest Fermilab data 
support this assumption. However, we need to further study the dependence 
of the hardness of fragmentations on both the beam atomic number. A, and 
the beam energy. 

Charmonium productions provide us with very interesting experimental 
data. The uncertainty from the fudge factors needs to be addressed in order 
to calculate the absolute yield of charmonium states. 

Another open problem that needs future investigation is the effect of 
the energy loss. It has been suggested that, due to the energy loss of charm 
quarks, the large initial open charm production could be kinematically shifted 
to the lower-mass region in dilepton spectra in ^4 -|- ^4 coUisions. To evalu- 
ate this possibility, parton cascade simulations will be needed. Such models 
are currently being formulated, and it will soon be possible to carry out de- 
tailed numerical simulations. The energy-loss effect could have interesting 
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consequences on the PHENIX search for direct thermal dileptons. In ad- 
dition, a detailed analysis on dilepton spectra from open charm decays in 
A + A collisions could provide information on the still-controversial energy- 
loss mechanism in the QGP. 
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Appendix 

In the appendix, we shall consider the impact-parameter-dependence of nu- 
clear shadowing effects. 



For the shadowing used in HIJING model W^, the density function is 



Ta/A{x,Q'^,r) = TA{r)fa/N{x,Q'^)SA{x,r) , (A.l) 

where fa/N{x,Q'^) is the parton structure function in a nucleon, r is the 
transverse vector of parton a inside nucleus A, TA{r) is the transverse density 
function of the nucleus, and for a uniform-sphere nucleus 



T^(f) = {3A/27tRI)^1 - ryR\ . (A.2) 

r) describes shadowing effects with the following parameterization, 
which is independent of the parton species and the scale: 

SA{x,r) = — 

= 1 + 1.19 In^/^ A[x^ - 1.5(xo + xl)x'^ + 3xoXlx] 



■ 1.08(AV3 - 1) ■ 

"^^^^ - ln(A + l) ^. 



(A.3) 



where 



xo = 0.1, XL = 0.7 , (A.4) 
4 



a^(f) = 0.1(^1/3 - l)-yi - rVi?i . (A.5) 
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For central p + Au collisions, we should take the average over r < ro( 
1.12/m), where ro <^ Ra{^ A^/'^tq). Therefore 



(A.6) 



and the integral over the transverse vector in Equation ( [4.28| ) gives the fol- 
lowing factor: 



d^rTp{r)TAu{r)Sa/p{x, r)SAu{x, r) 



2^ 



(A.7) 



Note that there is no shadowing on the proton projectile in p + Au collisions. 

The other shadowing functions from Eskola ||166|| depend on and are 
also different for valence quarks, sea quarks and gluons. The density function 



is: 



^a/A(a;,Q^r) = TAir)fa/N{x,Q'^) 

ATA{f^ 



X 1 



1 - Sa/A{X,Q^) 



If we average over the impact parameter, we get 

fa/A{x,Q'^) = J (frTa/A{x,Q'^,r) 

= Afa/N{X,Q')S,/A{X,Q') , (A.9) 

which is precisely the definition of the impact-parameter-independent shad- 
owing function Sa/AiXjQ"^)- For central p + Au collisions, the integral over 
the transverse vector yields 

ATaJ^ 



<frTp{r)TAu{^ll- 



27rrg 



4 



-^Sb/Au{.X, Q'^ )-- 



. 1- 



l — Sb/Auix, 



(A.IO) 
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From equations ( [A.7| ) and ( [A.10| ), the overall shadowing factor for central 
p + Au collisions is 

Sau{x,^) (A.ll) 

for HIJING's case, and 

[A/^S,/Au{x,Q^) -!/?>] (A.12) 

for Eskola's case. 
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